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ABSTRACT 

Observations of interstellar gas-phase and solid-state species in the 2.4-200 /im range obtained 
with the spectrometers on board the Infrared Space Observatory are reviewed. Lines and bands 
due to ices, polycyclic aromatic hydrocarbons, silicates and gas-phase atoms and molecules 
(in particular H 2 , CO, H 2 0, OH and CO2) are summarized and their diagnostic capabilities 
illustrated. The results are discussed in the context of the physical and chemical evolution of 
star-forming regions, including photon-dominated regions, shocks, protostellar envelopes and 
disks around young stars. 

Key Words: infrared: spectroscopy; star formation; interstellar molecules; interstellar dust; 
circumstellar matter 
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1. INTRODUCTION 



During the formation of stars deep inside molecular clouds, the surrounding gas and dust 
become part of the infalling envelope feeding the central object. The protostars themselves are 
extinguished by hundreds of magnitudes, so that the circumstellar gas and dust are the only 
tracers of the physical processes happening inside. A study of their characteristics is therefore 
key to understanding the origin of stars. Part of this gas and dust ends up in the rotating 
disks surrounding young stars, and forms the basic material from which icy planetesimals, 
and ultimately rocky and gaseous planets, are formed. Spectroscopic surveys of star-forming 
regions at different evolutionary stages therefore also provide quantitative information on the 
building blocks available during planet formation. The Infrared Space Observatory (ISO) 1 has 
provided the first opportunity to obtain complete infrared spectra from 2.4-200 fj,m unhindered 
by the Earth's atmosphere. 

At the temperatures characteristic of star-forming regions, most of the radiation emerges at 
mid- and far-infrared wavelengths, the majority of which is blocked from Earth. Most young 
stellar objects (YSOs) have been found through ground-based infrared surveys and the In- 
fraRed Astronomical Satellite (IRAS). The spectra of the coldest protostellar objects (ages of 
~ 10 4 yr since collapse began) peak around 100 fim and such sources are best studied with 
ground-based submillimeter telescopes and the longest wavelength instruments on ISO. Once 
the dense envelopes start to dissipate due to the effects of outflows, the objects become de- 
tectable at mid-infrared wavelengths, around ages of ~ 10 5 yr, and are accessible to the shorter 
wavelength instruments on ISO. The outflows also create shocks in the surrounding molecu- 
lar cloud, which emit brightly in mid-infrared lines. If ultraviolet (UV) radiation from the 
young stars can escape, it heats the neighboring gas and photodissociates molecules, creating 
so-called photon-dominated or photodissociation regions (PDRs). The tremendous range in 
physical conditions in these different phases, with densities from 10 4 to 10 13 cm -3 and tem- 
peratures from 10 to 10,000 K, are reflected in the abundances and excitation of the atoms and 
molecules. The ISO spectroscopic data therefore contain a wealth of information with which 
to unravel the physical and chemical processes during star- and planet formation. 

ISO's main contributions have been in the following areas. First, the complete wavelength 
coverage has provided an unbiased overview of the major gas- and solid-state species in star- 
forming regions, including an inventory of the reservoirs of the major elements (C, O, N, ...). 
Indeed, the identifications of many species have become much more secure because of the broad 
spectral range without gaps due to the Earth's atmosphere. Second, ISO was particularly well 
suited to study the physical structure of warm gas in the important 100-2000 K regime, which 
is difficult to probe with ground-based telescopes: (sub-)millimeter CO line observations trace 
cold (< 100 K) gas, whereas near-IR H2 and optical atomic lines refer to much hotter (> 2000 
K) gas. Third, ISO has quantified the energetics of star-forming regions, by measuring directly 
the photoelectric heating efficiency as well as the total cooling power in all relevant lines. 

ISO was a project of the European Space Agency with instruments funded by ESA member states (especially 
the PI countries: France, Germany, The Netherlands, and the United Kingdom) and with the participation of 
ISAS and NASA (Kessler et al. 1996, 2003). 



3 





6 


10 


15 














1 / 














3>i 








mm so 



140 iBH 150 



9 10 30 



Figure 1 (color). ISOCAM composite 7 and 15 /jm image of the p Oph molecular 
cloud. The bright filaments in the west trace the hot dust illuminated by the B2V star 
HD 147889 located just outside the map. The dense cloud containing IRAS 16293 
-2422 is to the east. The dark patches are very dense cores which are optically thick 
even at 15 /im, whereas the bright point sources throughout the image are low-mass 
YSOs. Most of the extended 7//m emission is due to PAHs (Abergel et al. 1996). 
Spectra of various objects in p Oph are indicated on the side. Left Top: PHOT-S 
spectrum of WL 6 (Gurtler et al., unpublished); Right Top: CAM-CVF spectrum of 
p Oph W (Boulanger et al. 2000); Left Bottom: LWS spectrum of IRAS 16293 -2422 
(Ceccarelli et al. 1998a); Right Bottom: SWS spectrum of Elias 29 (Boogert et al. 
2000b). 



Finally, the sensitivity of the instruments has allowed much larger samples to be probed, so 
that systematic trends and characteristics at each stage of evolution could be determined. 
Much of this progress would not have been possible without dedicated laboratory astrophysics 
and theoretical chemistry studies to provide the basic molecular and solid-state data needed 
to analyze the ISO spectra. 

This review summarizes the spectroscopic results from all four ISO instruments (see Table 1 
for characteristics). The bulk of the data come from the Short Wavelength Spectrometer 
(SWS) and Long Wavelength Spectrometer (LWS), but relevant spectra obtained with the 
Camera-Circular Variable Filter (CAM-CVF) and Photometer-Spectrometer (PHOT-S) are 
included. It is heavily biased toward ISO results only, but the science case builds on pioneering 
observations performed prior to ISO with ground-based and airborne infrared telescopes, in 
particular the Kuiper Airborne Observatory (KAO) (Haas et al. 1995), and with the IRAS 
and the InfraRed Telescope in Space (IRTS) missions. A historical review of those data in the 
context of the ISO results has been given by van Dishoeck & Tielens (2001). Also, more recent 
results from complementary ground-based facilities and new scientifically related space missions 
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Table 1. Spectrometers on the Infrared Space Observatory 



Instrument 


Wavelength 


Resolving power 


Aperture 


Ref. 




range (/xm) 


(A/AA) 


(arcsec) 




sws 


2.4-45.2 


1500 a 


14 x 20 - 20 x 33 


1 




11.4-44.5 


30000 


10 x 39 - 17 x 40 




LWS 


43-197 


200 


~ 80 6 


2 






10000 


~ 80 b 




CAM-CVF 


2.3-16.5 


35-50 


12 x 12 


3 


PHOT-S 


2.5-5, 5.8-11.6 


90 


24 x 24 


4 



a Varying with wavelength from 1000 — 2500; full spectral scans have usually been obtained 
with a factor 2 lower spectral resolution. 

b Diameter of circular aperture. 

References. 1. de Graauw et al. 1996a; 2. Clegg et al. 1996; 3. Cesarsky et al. 1996; 4. 
Lemke et al. 1996 

such as the Submillimeter Wave Astronomical Satellite (SWAS) and ODIN are mentioned only 
in passing. 

Due to space limitations, this review focuses on spectroscopy of galactic molecular clouds and 
the immediate surroundings of YSOs 2 . Figures 1-3 illustrate the striking variety of spectral 
features associated with active star-forming regions. ISO imaging has also contributed to 
our understanding of the interstellar medium, especially through spatial variations in grain 
components (e.g., Abergel et al. 2002), the discovery of infrared dark clouds (Bacmann et al. 
2000, Omont et al. 2003) and studies of the coldest pre-stellar cores (Ward-Thompson et al. 
2002, Krause et al. 2003), but details are not covered here. This paper also does not include 
ionized regions of the interstellar medium, nor the circumstellar material around evolved (post- 
) AGB stars and planetary nebulae (Molster & Waters 2003) . A general overview of interstellar 
dust has been given by Draine (2003) whereas extragalactic ISO spectra have been summarized 
by Genzel & Cesarsky (2000). Even with these limitations, it is not possible to refer to all 
ISO spectroscopy results. Many of them have been summarized in dedicated ISO conferences, 
including Heras et al. (1997), Yun & Liseau (1998), Waters et al. (1998), Cox & Kessler (1999), 
d'Hendecourt et al. (1999) and Salama et al. (2000). 



The term 'protostar' refers strictly only to objects which derive most of their luminosity from accretion and 
have L su bmm/Lbol > 0.5% and M env /M* > 1 (Andre et al. 2000). The term 'young stellar object' refers to 
all embedded stages of star formation and associated phenomena such as outflows. In this review, the term 
protostar will be used loosely to also refer to the more evolved embedded stages, especially in the context of 
high-mass objects. 
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Figure 2 (color). CAM-CVF image of the HH 1-2 region by Cernicharo et al. 
(2000b), illustrating a wide variety of types of spectra, ranging from ice absorptions 
toward the class protostar VLAl (spectrum b2) to silicate emission from a disk 
around the Cohen-Schwarz star (spectrum b3) and gas emission from shocks (spectra 
bl, b4, b5). Both low-velocity shocks (v ~ 10 — 15) km s _1 as traced by the strong 
H2 lines in panel bl and high- velocity shocks (v ~ 100 — 140 km s _1 , as traced by 
the [Ne II] and [Ne III] lines in panel b4 are seen. Note the complete absence of any 
continuum emission due to hot dust and PAH emission in both panels. The color 
image is a [Si II] 6717 A optical image, whereas the black contours indicate the 7 //m 
CAM image. The intensity scale is in mJy 

2. INFRARED SPECTROSCOPIC FEATURES 

The mid- and far-infrared wavelength ranges are spectacularly rich in spectral features, many 
of which probe unique species which cannot be observed at any other wavelength. The wealth 
of sharp gas-phase lines and broad solid-state bands is illustrated in the SWS spectrum of one 
of the brightest mid-infrared sources in the sky, the Orion-KL region (Figure 3). Each of these 
features, summarized in Table 2, has its own diagnostic capability. 

Atomic lines: Fine-structure transitions within the lowest electronic term of most astrophys- 
ically relevant atoms and ions occur at mid-infrared wavelengths. Many of them are found in 
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Table 2. Selected mid- and far-infrared spectral features observed by ISO 



Category 


A 


Species / 
Line 


Diagnostic 0, 


PDR Shock 


Embedded 
YSO 


Disk 


Atoms 


25.2 


[S I] 


shock- vs. 


+ 


- 


- 




34.8 


[Si II] 


photon heating 


+ + 


- 


- 




63.2 


[O IJ 




+ + 


+ 






145.5 


fO 11 

L J 




+ + 


+ 






157.7 


[C II] 




+ + 


+ 


? 


H 2 


6.9 


7^5 


mass and temperature 


4- 4- 








8.0 


6^4 


of warm gas, 


4- 4- 








9.7 


5^3 


shock- vs. 


4- 4- 


+ 






12.2 


4^2 


photon heating 


4- 4- 


4- 






17.0 


3->l 




4- 4- 


4- 


? 




28.2 


2^0 




-f -f 


+ 


? 


HD 


19.4 


6^5 


[D]/[H] 


+ 


- 


- 




112.0 


l-»0 




+ 


- 


- 


Gas-phase 


6.0 


H 2 


temperature 4- density, 


4- 


4- 




Molecules 


7.7 


CH 4 


ice evaporation, 




-|- 






13.7 


C2 H2 


organic chemistry, 




4- 






14.0 


HCN 


depletion 




4- 






15.0 


C0 2 




+ 


+ 






104.4 


CO 25^24 




+ 


+ 


- 




108.1 


o-H 2 2 2 i - lio 




+ 


+ 


- 




119.3 


oh 2 n 3/2 2^2 




+ 


+ 


- 




1 Q A A 
13U.4 


(~^r~\ on , 1 




+ 


+ 






138.5 


d-HqO Si** — 2o9 

f 1 1 Z u lo ^oz 




-f 








162.8 


CO 16^15 




+ 


+ 


_ 




174.6 


0-H2O 3 03 - 2i2 




+ 


+ 






179.5 


O-H2O 2l2 - 1 1 




+ 


+ 


- 




186.0 


CO 14^13 




+ 


+ 


- 


PAHs 


3.3, 6.2, 7.7, 


C-H and 


carbonaceous material, 


4- 




+ 




8.6, 11.3, 12.7, 


C-C modes 


UV radiation 


+ 


_ 


+ 




14.2, 16.2 






+ 


- 


+ 


Silicates 


9.7 




bulk of dust 




-4- 


4- 


(Amorphous) 


18.0 








4- 


-f- 


Silicates 


1 1 Q 1 a A OQ o 

11.3, 16.4, 


Forsterite 


mineralogy, 






+ 


(Crystalline) b 


27.7, 33.8, 69 




grain growth and 






4- 




18.5, 21.5, 24.5 


Enstatite 


processing/heating, 






+ 




8.6 


Silica 


solar system 






4- 




65 


Diopside 


connection 




? 


? 


Oxides 


11.6 


AI2O3 


solar system 


- 


- 


+ 




23 


FeO 


connection 


- 


- 


+ 


Sulfides 


23.5 


FeS 




- 




+ 


Carbonates 


92.6 


Calcite 






? 




Ices 


4.27 


CO2 


volatile solids, 




+ 






4.38 


13 co 2 


organic chemistry, 




+ 






4.67 


CO 


thermal history, 




+ 






6.0, 13 


H 2 


solar system 




+ 






6.85 


CH 3 OH + NH+? 


connection 




+ 






7.7 


CH 4 






+ 






15.2 


C0 2 






+ 






44, 63 


H 2 (cryst.) 






+ 


+ 



a Diagnostic properties of category of species; individual species or lines probe a subset of these properties. 



Position may vary depending on composition 




H 2 0-ice 



0.1 I I I I I I I I I I I l_ 

2 3 4 5 10 20 30 40 

wavelength [fj,m] 



Figure 3. SWS grating scan of the Orion 'Peak 1' shock, showing a rich forest of H 2 
lines and other features (Rosenthal et al. 2000). 



H II regions where they are used to probe the hardness of the radiation field, but some lines 
are observed in the neutral clouds discussed here, e.g., [S I] 25.2, [Si II] 34.8, [O I] 63.2 and 
145.5 and [C II] 157.7 /iin 3 . These are key diagnostics of photon- versus shock-heating. H I 
recombination lines are often seen in mid-infrared spectra of ionized gas, but not in those of 
neutral clouds. 

H2 and HD pure rotational lines: The lowest transitions of the most abundant molecule in 
the universe, H 2 , occur at mid-infrared wavelengths, and probe directly the bulk of the warm 
gas and its temperature. Most of the deuterium in dense clouds is in HD, whose fundamental 
transition lies at 112 fim. 

Gas-phase bands: Fundamental vibrational transitions of important molecules such as H 2 0, 
CH 4 , C2H2, HCN and CO2 occur at mid-infrared wavelengths. This is the only way to observe 
symmetric molecules like CH4 and C2H2 which have no dipole moment and thus cannot be 
observed through rotational transitions at millimeter wavelengths. Also, CO2 and H2O are 
so abundant in the Earth's atmosphere that their interstellar lines can only be detected from 

The following nomenclature is used: X refers to the atom or molecule as a chemical species, [X] indicates the 
abundance of the element in all forms with respect to total hydrogen, and [X I] 100 /im denotes its forbidden 
mid- or far-infrared line. 
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Figure 4. Normalized SWS spectra toward massive protostars showing absorption 
by various species. Left: CO2 v-i and CO v\ bands; Right: H2O ^2 band. The spectra 
in the right panel have been shifted vertically for clarity. Gas-phase molecules show 
sharp unresolved ro-vibrational lines in a characteristic pattern, whereas ices have 
a single broad absorption (Based on van Dishoeck 1998, Helmich et al. 1996, and 
Boonman & van Dishoeck 2003). 

space. The pure rotational transitions of H2O, OH and CO lie at far-infrared wavelengths. 
Together, the vibrational and rotational lines are good probes of the physical conditions of the 
gas, its cooling and its chemistry. 

PAHs: The C-C and C-H stretching and bending modes of polycyclic aromatic hydrocarbons 
(PAHs) at 6.2, 7.7, 8.6, 11.3, ... \xm dominate the mid-infrared spectra of many objects and 
are indicators of the presence of complex carbonaceous material excited by UV radiation. 

Solid-state vibrational bands: The characteristic vibrational bands of ices, silicates, oxides, 
carbides, carbonates and sulfides occur uniquely at infrared wavelengths. They can be used as 
diagnostics of heating as well as grain growth. Solid-state species can be distinguished from 
gas-phase molecules because their bands lack the characteristic ro-vibrational structure and 
are broadened (see Figure 4). 

Not all features listed above have been observed with all ISO instruments. In general, the 
detection of intrinsically narrow gas-phase lines in emission or absorption against a strong 
continuum such as Orion-KL required the highest spectral resolution provided by the Fabry- 
Perots of the SWS and LWS. For most other sources, sensitivity considerations dictated the 
use of the grating. The SWS spectral resolving power of 1000 — 2500 was particularly well 
matched to that needed to resolve the solid-state features and PAH bands, and analyse their 
substructure. At the low spectral resolving power provided by CAM-CVF and PHOT-S, the 
PAHs and broad solid state bands are readily detected in much weaker sources than accessible 
with the SWS, but the analysis is limited because the profiles are unresolved. Strong gas-phase 
lines are only detected with CAM-CVF if the continuum emission is nearly absent, such as in 
shocks or PDRs. 
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The main limitations of the ISO data are their relatively low spatial and spectral resolu- 
tion compared with modern ground-based telescopes. ISO's telescope was only 60 cm, with 
spectrometer apertures ranging from 20" to > V . This means that the different physical 
components associated with star formation (envelopes, outflows, disks) are often unresolved 
and blurred together. For extended sources, aperture changes cause jumps in the fluxes with 
wavelength. The spectral resolution was generally insufficient to obtain kinematic information. 

In Sections 3-7, each of the species from Table 2 is discussed in detail, focussing on their 
spectroscopy. Subsequently, the features are put in the context of the different physical re- 
gions associated with star formation — PDRs, shocks, protostellar envelopes and circumstellar 
disks — illustrating their diagnostic capability and ISO's contribution to our understanding of 
physical and chemical processes. 

3. GAS-PHASE MOLECULES AND THEIR CHEMISTRY 

Although hampered by limited spectral resolution, the SWS and LWS have observed many 
gas-phase species. ISO's main contributions have been of three types: (i) Detections of new 
interstellar molecules not seen prior to ISO (e.g., CO2, HF, CH 3 ; see Figure 5); (ii) Infrared 
detections of molecules seen previously at other wavelengths (e.g., H2O, SO2, HD); (hi) Infrared 
bands of molecules observed previously from the ground or KAO (e.g., H 2 , CH 4 , C3, C2H2, 
OH). 

3.1. New Detections 

Searches for gaseous CO2 formed one of the main astrochemical goals of the SWS. This molecule 
is potentially one of the more abundant carbon- and oxygen-containing species, but ground- 
based searches were limited to the chemically related species HOCO + . Detection of the CO2 ^2 
band at 14.98 /xm was reported by van Dishoeck et al. (1996) and confirmed by van Dishoeck 
(1998), Dartois et al. (1998b) and Boonman et al. (2003b, c) in a variety of star-forming regions. 
The inferred abundances with respect to H 2 are typically ~ 2 x 10~ 7 , up to two orders of 
magnitude lower than those of solid CO2 in the same regions (see §4). The line profiles 
indicate that the abundances can be a factor of 10 higher in warm gas at T > 300 K, whereas 
it is at least an order of magnitude lower in shocks. Its chemistry is not yet fully understood. 

Gas-phase CH 3 has been detected with the SWS through its bending mode lines at 16.0 and 
16.5 /im along the line of sight toward Sgr A* with an abundance of ~ 10 -8 (Feuchtgruber et 
al. 2000). This molecule is an important building block in the formation of small hydrocarbon 
molecules. Together with observations of CH, CH4 and C2H2 along the same line of sight, it 
forms a good test of basic gas-phase astrochemistry networks in diffuse and translucent clouds. 

HF is interesting because it is the main reservoir of fluorine in dense molecular clouds and 
thus a direct measure of the fluorine depletion, supposedly characteristic of other species. The 
detection of the HF J=2-l line at 121.7 fim in absorption toward Sgr B2 by Neufeld et al. 
(1997) implies an abundance of ~ 3 x 10~ 10 , indicating that only ~ 2% of the fluorine is in 
the gas-phase. 



10 



■ fH 

in 

a 



IE 
03 



i> 

HI— 1 
-Hi 

OS 




0.9 - 



0.8 



1 1 1 


co 2 






: 1 

i 


i 



1.03 - 



14. Q 15.2 
Wavelength (/^m) 




200 
V LSR (km/a) 



-200 

V 1SR (km/a) 



200 



Figure 5. New gas-phase molecules detected by ISO. The light line in each panel 
indicates the model spectrum. Top Left: CH3 U2 Q-branch toward Sgr A* (Feucht- 
gruber et al. 2000); Top Right: CO2 vi band toward AGFL 2136 (van Dishoeck et al. 
1996, Boonman et al. 2003b); Bottom Left. HF J = 2 - 1 toward Sgr B2 (Neufeld 
et al. 1997); Bottom Right: HD J = 1 — line toward the Orion Bar (Wright et al. 
1999). 

A new band at 57.5 /xm has been reported by Cernicharo et al. (2002) toward Sgr B2 and may 
be due to either C4 or C4H. If ascribed to C4, this would be the first detection of this molecule 
in outer space, at an abundance comparable to that of C3. The latter molecule has been seen 
toward Sgr B2 through many LWS lines by Cernicharo et al. (2000a), with an abundance of 
~ 3 x 10" 8 . 

Several other new molecules have been detected in the envelopes around late-type stars and 
proto-planetary nebulae. Although not formally part of this review, the exciting SWS discov- 
eries of CeHg (benzene), CgE^ and C4H2 (Cernicharo et al. 2001a) and C2H4 (Cernicharo et 
al. 2001b) toward CRL 618 should be mentioned. 

3.2. Other Molecules 



Most interstellar molecules with dipole moments are detected through their pure rotational 
emission lines at millimeter wavelengths. Complementary infrared absorption data provide 
valuable constraints on the physical structure and geometry of the region. Using ISO, this 
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Figure 6. Cartoon of the environment of young stellar objects and the origin of 
the various spectral features observed by ISO. The top diagram indicates the typical 
temperature structure of a protostellar envelope, with the critical temperatures for 
freeze-out, ice evaporation, and high-temperature chemistry indicated. 

has been possible for a few species, most notably HCN (Lahuis & van Dishoeck 2000; see also 
ground-based observations by Evans et al. 1991) and SO2 (Keane et al. 2001a). In both cases, 
the inferred abundances from the infrared data are two orders of magnitude higher than those 
obtained from the millimeter lines. Since the pencil-beam absorption data are more sensitive 
to the warm gas close to the protostar, the natural conclusion is that the abundances 'jump' 
by a factor of ~ 100 in the inner region of the envelope due to a combination of ice evaporation 
and high-temperature gas-phase chemistry (Boonman et al. 2001, see Figure 6). 

Molecules without a permanent dipole moment observed with ISO include CH4 (e.g., Boogert 
et al. 1998) and C2H2 (e.g., Lahuis & van Dishoeck 2000, Boonman et al. 2003c). Both species 
have been detected from the ground, but ISO allowed surveys in a larger number of sources. 
Analysis of the strong C2H2 Q— branch profiles gives rotational temperatures up to 1000 K. 
These temperatures correlate well with those of CO and other temperature tracers, indicating 
that C2H2 is a good probe of the inner envelope. 

Higher rotational lines of NH, NH2 and H3O" 1 " have been seen in absorption with the LWS 
toward Sgr B2 (Cernicharo et al. 2000a, Goicoechea & Cernicharo 2001a, b). Together with 
the detection of a wealth of far-infrared lines of hot NH3 (Ceccarelli et al. 2002a), this allows 
tests of the basic nitrogen gas-phase chemistry. The interpretation is complicated, however, by 
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the presence of many absorption components along the line of sight, including diffuse clouds, 
a warm dense envelope and a newly-revealed layer of hot shocked gas. Several other molecules 
including SH, H 2 + and CH2 have been searched for but not detected (Cernicharo et al. 2000a). 
The analysis of the full Sgr B2 Fabry- Perot spectrum is still in progress and may lead to further 
discoveries in this spectral 'goldmine. Higher-lying pure rotational transitions of CO + have 
been seen with the LWS toward the low-mass protostar IRAS 16293 -2422 by Ceccarelli et al. 
(1998b), where its high abundance indicates either the presence of a dissociative shock or the 
importance of UV or X-rays in the chemistry. 

3.3. H 2 and HD 

Interstellar H2 is commonly observed through UV absorption lines in diffuse gas, and through 
near-infrared emission lines in warm clouds, but only few mid-infrared pure rotational lines 
had been seen prior to ISO (e.g., Parmar et al. 1991). As discussed in detail in §7 and 8, the 
SWS has routinely observed many H 2 pure rotational lines in a variety of regions, providing an 
excellent new probe of their physical conditions. The detection of the fundamental J=1-^> 
line at 28.22 /urn was first reported by Valentijn et al. (1996) in the galaxy NGC 6946, and was 
subsequently seen in several galactic sources. This line probes the bulk of the warm gas at a 
temperature that is often lower than that inferred from the higher-lying H 2 lines. 

Like H 2 , interstellar HD has been observed through ultraviolet absorption lines prior to ISO, 
but only in diffuse clouds where the bulk of the deuterium is still in atomic form. The LWS has 
allowed the first detection of the fundamental HD J=l-0 line at 112 fim in emission in a dense 
cloud in Orion where most of the deuterium is in molecular form (Wright et al. 1999, see Figure 
5). In addition, the higher-lying HD J=6-5 line at 19.4 /um has been seen with the SWS in the 
Orion shock (Bertoldi et al. 1999). Together with observations of the pure-rotational H 2 lines, 
an accurate HD/H 2 abundance — and thus [D]/[H] abundance — can be derived. For both 
sources, a [D]/[H] ratio of (0.5 — 1.0) x 10 -5 has been found, up to a factor of two lower than 
the [D]/[H] ratio of (1.5 ±0.1) x 10" 5 measured in local diffuse clouds (Moos et al. 2002). The 
lower [D]/[H] ratio in an active star-forming region like Orion may be evidence of deuterium 
destruction by nuclear burning since its creation in the Big Bang. The HD 112 //m line has also 
been reported in emission in Sgr B2 (Polehampton et al. 2002) and in absorption toward W 
49 (Caux et al. 2002). In these latter cases, H 2 is not observed directly and the determination 
of [D]/[H] is complicated by the fact that alternative tracers of H 2 , in particular CO, may be 
significantly frozen out onto grains. 

3.4. H 2 0, OH and O 

Various H 2 lines have been seen prior to ISO from the ground and KAO at radio and mil- 
limeter wavelengths, but in virtually all cases these observations concern maser transitions. 
ISO has opened up a flood of data on thermal H 2 lines, either in emission or absorption, 
allowing much more accurate determinations of the H 2 abundance and excitation. Pure 
rotational emission and absorption lines at 25-200 fim have been seen with the LWS (e.g., 
Liseau et al. 1996, Cernicharo et al. 1997, Harwit et al. 1998) and SWS (Wright et al. 2000), 
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and vibration-rotation lines at 6 with the SWS (e.g., Helmich et al. 1996, Dartois et al. 
1998b, Gonzalez- Alfonso et al. 1998, 2002, Moneti & Cernicharo 2000, see Figure 4). Large 
variations in the H 2 abundances are found, ranging from < 10~ 8 in the coldest clouds to 
> 10~ 4 in warm gas and shocks, illustrating the extreme sensitivity of this molecule to the 
physical conditions. These variations can be explained by freeze-out of water at the lowest tem- 
peratures, ice evaporation at warmer temperatures ( ^ 90 K) and gas-phase reactions driving 
most of the oxygen into water at high temperatures (> 230 K) (see Figure 6, §9.2). 

The chemically related OH molecule has been detected through its pure rotational far-infrared 
transitions in various objects (see §8), most prominently in absorption toward Sgr B2 (Goi- 
coechea & Cernicharo 2002). Along this line of sight, [O I] 63 fim absorption has been seen 
as well, both with the LWS grating (Baluteau et al. 1997) and the LWS Fabry-Perot (Lis et 
al. 2001, Vastel et al. 2002). The conclusion from both sets of data is that the bulk of the 
gas-phase oxygen is in atomic form with 0/CO~ 3 — 9, consistent with previous KAO results 
for DR 21 by Poglitsch et al. (1996). Even higher O/CO ratios have been inferred for W 49N 
by Vastel et al. (2000) and for the cold cloud L1689N by Caux et al. (1999). Although care 
should be taken in analyzing highly optically thick, potentially self-absorbed [O I] emission 
lines, the inferred high atomic oxygen abundances are consistent with non-detection of O2 by 
the SWAS satellite (Goldsmith et al. 2000). 

4. INTERSTELLAR ICES 

Interstellar H2O ice was detected at 3 /xm by Gillett & Forrest (1973), and the study of ices has 
subsequently been pursued from the ground and with the KAO (e.g., Willner et al. 1982). ISO 
has doubled the number of detected ice bands, bringing the current total to nearly 40 features. 
Most of the bands are observed in absorption arising in the cold outer envelopes of high-mass 
protostars (see Figures 6 and 7), but ISO has also allowed the first glimpse of ices toward 
low-mass objects (Cernicharo et al. 2000b, see Figure 2). The absorptions are deep: the total 
integrated optical depth is often larger than that of the 9.7 and 18 //m silicate absorptions and 
can amount to more than 50% of the total integrated flux. Only a few 'windows' in between 
the ice bands remain to probe close to the protostar itself. Ice bands are also seen with the 
LWS, in absorption toward protostars at 44 /im (Dartois et al. 1998a) and in emission in 
disks around young stars at 44 and 63 /an (Malfait et al. 1998, 1999, Creech-Eakman et al. 
2002) and one shock (Molinari et al. 1999). In the latter cases, water ice is in its crystalline 
form rather than in the amorphous phase which is thought to dominate most of interstellar 
ice. In terms of abundances, the amount of ice (mostly H 2 0) can be comparable to that of 
the most abundant gas-phase molecule containing heavy atoms (CO), making it the second 
most abundant species after H2 in cold clouds. Thus, a good understanding of the processes 
involving ices is highly relevant to understanding the physical and chemical characteristics of 
star-forming regions. Conversely, because ices are such a major component and show large 
variations in abundances and profiles, they are particularly powerful diagnostics of changes 
in environment. Recent reviews on interstellar ices and summaries of ISO data are given by 
Schutte (1999), Ehrenfreund & Schutte (2000), Boogert & Ehrenfreund (2004), and Gibb et 
al. (2004). 
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Figure 7. SWS spectrum of the embedded massive protostar W 33A, illustrating 
the wealth of ice features observed by ISO (Gibb et al. 2000). Blow-ups of the 5-8 
/tm region are shown for W 33A, NGC 7538 IRS9 (shifted in flux for clarity) and Mon 
R2 IRS3. Note the shift in the 6.85 /im band toward longer wavelengths for warmer 
sources like Mon R2 IRS3 (Keane et al. 2001b). 

ISO's main contribution to the study of ices has been in two areas: (i) The first complete inven- 
tory of the main ice features; and (ii) Clear evidence that changes in profiles and abundances 
can trace the gradual heating of protostellar envelopes. 

4-1. Inventory of Ice Features 

Table 2 and Figure 7 show a number of features with their identifications (see Boogert & 
Ehrenfreund 2004 for a complete list). Several of them are securely identified on the basis 
of multiple vibrational modes and isotopic bands, such as H 2 0, CO, CO2 and CH3OH. For 
example, Dartois et al. (1998a) are able to fit five H 2 bands simultaneously. Cleanly isolated 
single bands such as that due to solid CH4 at 7.67 /im and OCS at 4.92 /jm can also be 
confidently assigned if they are well matched by laboratory spectra. The identification of 
other species like H 2 CO and NH 3 is more problematic because all their strong vibrational 
modes are blended with those of other species. 

High quality SWS profiles of the 6.0 and 6.8 /im bands — prominent toward all ice sources 
- show that they are clearly composites of several features (Schutte et al. 1996, Keane et al. 
2001b). A good fit to the 6.0 /im band requires not only H2O (known to be present from the 
3 /im band), but also H2CO at 5.85 /im and/or HCOOH at 5.83 /im to fit the blue wing and 
NH 3 at 6.15 /mi or some carbonaceous material to fit the red wing. The 6.85 /im feature is still 
not firmly identified, although NH^" is a plausible candidate and CH3OH a minor contributor 
(Demyk et al. 1998, Keane et al. 2001b, Schutte & Khanna 2003). An empirical decomposition 
of the profile shows two components with fixed position and widths, whose relative intensities 
change with temperature of the source (see Figure 7). 

The shape, exact position and intrinsic strength of each band depend on the chemical envi- 
ronment of the molecule. For example, the position of the CO2 stretching mode shifts by 20 
cm -1 and its width broadens by a factor 2 when CO2 is mixed with H2O ice rather than in 
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its pure form. Analysis of the band profiles has led to a picture of layered ices, in which some 
molecules (e.g., CO) are not mixed with H2O (Tielens et al. 1991). Ice profiles of crystalline 
material can be distinguished from their amorphous counterparts by being sharper and red- 
shifted. Laboratory simulations of astrophysically relevant ice mixtures have been essential 
to provide the basic data to quantitatively analyze the ice bands (e.g., Hudgins et al. 1993, 
Gerakines et al. 1995, Ehrenfreund et al. 1996, 1997a). 

A number of weak minor features have been found in the 7-8 /jm region. Schutte et al. (1999) 
and Keane et al. (2001b) suggest HCOOH and HCOO - or CH3HCO as plausible identifications. 
The presence of ions such as NH+, OCN" and HCOO" has been subject to considerable debate 
over the last 20 years, with the major controversy centering around the charge balance in the 
ice (not all counter ions have been found) and the constancy of the position of the observed 
ionic features. Schutte & Khanna (2003) show that heavier counterions have only weak spectral 
features and likely remain undetected, whereas the observed 6.85 /im shift is well reproduced 
in the laboratory by heating NH^-containing ices. 

Upper limits can be just as interesting as detections. The region between 5.0 and 5.8 /im 
(Figure 7) is remarkably void of features and places limits on minor ice constituents. For 
example, the absence of the 5.62 /im band of the simplest amino acid glycine puts a limit on 
its abundance of less than 0.3% with respect to H2O ice (or < 3 x 10 -7 with respect to H2). 
Boudin et al. (1998) have put significant upper limits on C2H6, C2H5OH and H^O^- Solid 
HNCO, the likely precursor of OCN - and widely observed in the gas in hot cores, absorbs at 
4.42 [xva. and is undetected down to 0.5%. In these regions, the 'confusion' limit has certainly 
not yet been reached and higher S/N data may reveal new weak bands. Similarly, very few 
bands have been seen superposed on the strong continuum longward of the 9.7 fim silicate 
band, with the strong CO2 bending mode at 15 /im being the notable exception. 

Searches for absorption by more complex solid organic material are of obvious interest as a 
link with the large organic gas-phase molecules observed toward high-mass protostars. PAHs 
have been proposed to contribute to the red wing of the 6.0 fim absorption. Refractory organic 
residues, put forward by Greenberg et al. (1995) as resulting from UV processing of ices, have 
been suggested as an alternative explanation for the 6.0 ^m excess absorption (Gibb & Whittet 
2002). 

Two very important molecules, O2 and N2, are largely "invisible" because their vibrational 
modes are dipole forbidden. When surrounded by other molecules, however, the transitions 
become weakly allowed. Searches for the N 2 band at 4.295 /im are hampered by overlap with 
the strong solid CO2 stretching mode (Sandford et al. 2001). Upper limits on solid O2 feature 
at 6.45 ^m, combined with an analysis of the effects of O2 on the solid CO profiles, indicate 
that solid O2 contains less than 6% of the total interstellar oxygen budget (Vandenbussche et 
al. 1999). 

The isotopic species 13 C02 and 13 CO have been detected and provide an independent deter- 
mination of the [ 12 C]/[ 13 C] abundance gradient as a function of galactocentric radius as well as 
constraints on grain shape effects which complicate the analysis of the main isotopes (Boogert 
et al. 2000a). The potential detection of solid HDO by Teixeira et al. (1999) has not been 
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confirmed by a re-analysis of the ISO data nor by independent ground-based data (Dartois et 
al. 2003). 

4-2. Ice Abundances and Chemistry 

H 2 ice is clearly dominant, with abundances of 10 -5 — 10 -4 with respect to H2, thus containing 
a significant fraction of the oxygen (van Dishoeck 1998). The combined contribution of CO 
and CO2 is typically 10-30%, but may be as high as 60%. For some sources, especially the 
most massive protostars, other carbon-bearing species (CH3OH, HCOOH, CH 4 and H 2 CO) 
may add an additional 20-30%, bringing the total abundance up to ~50% of that of H2O 
(Gibb et al. 2004). Nitrogen appears much less abundant in the ice. Although the amount 
of NH 3 ice is under debate for some sources, it is certainly less than 10% and often less than 
5% of that of H 2 ice (e.g., Taban et al. 2003). The NHj abundance may be of order 10% 
(Schutte & Khanna 2003), but this still makes the identified nitrogen content at least a factor 
of 4 less than that of carbon. 

A striking aspect of the ice composition is its simplicity. Although there is some observational 
bias against detection of low-abundance complex molecules in ices, all of the identified species 
are those expected from hydrogenation and oxidation of the main species to be accreted from 
the gas, i.e., O, C, N and CO (Tielens & Hagen 1982). Various groups (e.g., Charnley & 
Tielens 1987, Stantcheva et al. 2002) have attempted to model this grain surface chemistry 
quantitatively with varying degrees of success. 

The SWS has provided detailed inventories of ices in about 30 sources. The majority of them 
are high-mass protostars, with only a limited number of intermediate mass YSOs (LkHa 225: 
van den Ancker et al. 2000c; Elias 18, HH 100, R CrA: Nummelin et al.2001; AFGL 490: 
Schreyer et al. 2002) and low-mass sources (Elias 29: Boogert et al. 2000b; L1551 IRS5: White 
et al. 2000; T Tauri: van den Ancker et al. 1999; see also van den Ancker 2000 for unpublished 
SWS results). Using the higher sensitivity provided by CAM-CVF, Alexander et al. (2003) 
surveyed 42 low luminosity sources (L <100 L Q ), whereas PHOT-S observed a handful of cases 
(e.g., Feldt et al. 1998, Gurtler et al. 1996, 1999, 2002). Although the lower spectral resolution 
does not allow a decomposition of the intrinsically blended 6.0 and 6.8 /im profiles, their overall 
spectra look remarkably similar to those of the higher luminosity counterparts. The 6.8 /im 
feature clearly correlates with H 2 0, confirming that it is largely due to an ice. Significant 
variations are found between different clouds. For example, YSOs in the Chamaeleon I cloud 
show very little ice absorption overall, whereas those in p Oph appear to have a lower H 2 ice 
abundance and a larger scatter in the C0 2 /H 2 abundance ratio compared with, for example, 
Serpens and CrA. p Oph may be peculiar because it has a large foreground layer of molecular 
gas with Ay > 10 mag which may be relatively poor in ices due to its lower density and larger 
exposure to UV radiation (Boogert et al. 2002). 

4-3. Heating and Processing of Ices 

The 15 /im bending mode of solid C0 2 toward high-mass protostars shows considerable vari- 
ation from source to source, with a pronounced double-peak structure and a red wing seen in 
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some sources (de Graauw et al. 1996b, Gerakines et al. 1999). Extensive laboratory simulations 
show that these variations can be well reproduced by gradual heating of a mixture of H2O, 
CO2 and CH3OH ices in roughly equal proportions (Ehrenfreund et al. 1998, Dartois et al. 
1999). At higher temperatures, the laboratory ices appear to segregate, leading to near-pure 
H 2 and CO2 structures. Although this explanation seemed far-fetched at first, it is corrob- 
orated by a wealth of other evidence. In particular, the profiles of solid 13 C02 at 4.38 /im 
(Boogert et al. 2000a), H2O at 3 /im (Smith et al. 1989), and the 6.8 /jm feature (Keane et al. 
2001b) show similar changes due to heating. Moreover, these changes are well correlated with 
increased gas/solid ratios derived from ISO (van Dishoeck & Helmich 1996, Boonman & van 
Dishoeck 2003, see §9.2) and submillimeter data (van der Tak et al. 2000a) as well as increased 
dust color temperatures derived from the 45/100 /jm ISO flux ratios (Boogert et al. 2000a). 
Indeed, the predictive power of these heating indicators is large: if one of them is found, so 
are the others. The emission data refer to the entire source rather than absorption along the 
line of sight, indicating that the entire envelope undergoes 'global warming' rather than lo- 
cally. While geometry (e.g., where does the mid-infrared continuum become optically thick?) 
and line-of-sight effects (e.g., absorption looking down the outflow cone) may play a role in 
individual cases, they cannot be responsible for these general trends. As discussed further in 
§9.2, the most likely scenario is that the higher temperatures are due to a decreased envelope 
mass relative to the source's luminosity, perhaps caused by the dispersion of the envelope with 
time. 

The presence of some features, e.g., that due to OCN~ at 4.62 /im, has often been cited as 
an indicator of 'energetic' processing, i.e., changes induced in the ice matrix composition or 
profiles due to UV or particle irradiation rather than thermal heating (e.g., Pendleton et al. 
1999). There have been extensive laboratory studies of these processes (see Gerakines et al. 
1996 for UV; Moore & Hudson 1998 for energetic particle bombardment), but no 'smoking gun' 
has yet been found observationally. In particular, van Broekhuizen et al. (2004) show that the 
OCN - data can be quantitatively reproduced by thermal acid-base reactions of HNCO with 
NH 3 starting at temperatures as low as 10 K. 

5. POLYCYCLIC AROMATIC HYDROCARBONS 

Broad, strong infrared emission features at 3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 /jm from a variety 
of astronomical sources have been known since the 1970's, and are commonly known as the 
Unidentified InfraRed (UIR) bands. These features coincide with the vibrational modes of 
aromatic materials (Duley & Williams 1981), with the preferred identification that of free- 
flying Polycyclic Aromatic Hydrocarbons (PAHs) (Leger & Puget 1984, Allamandola et al. 
1985). The fact that the bands are detected far from any illuminating source with roughly the 
same intensity ratios as close to bright stars proves that the carriers are excited non-thermally 
by single UV photons (Sellgren 1984). Indeed, both CAM-CVF, PHOT-S and the IRTS have 
shown that the mid-infrared spectra (and thus also the IRAS 12 fim emission) of even the 
most diffuse interstellar medium are dominated by PAH features (Mattila et al. 1996, Onaka 
et al. 1996, Boulanger et al. 2000, Kahanpaa et al. 2003). Any large (>100 A) grain does not 
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Figure 8 (color). SWS spectrum of the Orion Bar PDR, illustrating the richness 
of PAH features. Various atomic and molecular lines are seen as well, originating in 
the PDR and in the Orion H II region along the line of sight. The spectrum is taken 
at the position of the peak of the H2 v = 1 — S(l) emission (Peeters et al. 2004b). 



become hot enough under these conditions to emit at mid-infrared wavelengths, so the carriers 
must be small, of order 10 A, corresponding to PAHs with 30 to a few hundred carbon atoms 
containing up to 10-15% of the available interstellar carbon. Larger PAH clusters are seen 
through the plateaus underlying the discrete features. PAHs have also been seen in absorption 
at 6.2 (im in the diffuse medium toward objects in the Galactic Center and several dusty 
late- type Wolf-Rayet stars (Schutte et al. 1998). 

The main contributions of ISO in this area have been to (i) Allow a systematic analysis of the 
spectral characteristics of PAHs, revealing many new weak features, subbands and shoulders; 
(ii) Survey the similarities and differences between the features for a wide variety of sources, 
relating them to physical and chemical changes; and (iii) Extend PAH studies to more diffuse 
regions exposed to much lower intensity radiation fields (see above). PAHs have been detected 
in many interstellar ISO spectra, indicating that they are ubiquitous, abundant and able to 
survive in a wide range of conditions. Deeply-embedded YSOs and shocks are the only types 
of regions where PAH emission is not seen. Detailed reviews of ISO PAH results have been 
given by Tielens et al. (2000) and Peeters et al. (2004b). 
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5.1. PAH Spectroscopy 



Figure 8 shows the rich PAH spectrum of the Orion Bar PDR. In addition to the main features, 
weaker bands or substructures can be found at 5.2, 5.7, 6.0, 7.2-7.4, 7.6, 7.8, 8.2, 10.8, 11.0, 
11.2, 12.0, 12.7, 13.2, 14.5 and/or 16.4 ^m, with a weak plateau between 15 and 20 /jm (see 
also Moutou et al. 2000, van Kerckhoven et al. 2000). Assignments in terms of C-H and C-C 
stretches are indicated. The C-H out-of-plane bending modes at 11-14 //m are particularly 
useful to characterize the PAH edge structure since their position depends on the number of 
adjacent H atoms bonded to neighboring C atoms on the ring. Solo-CH groups at 11.3 fim 
are indicative of straight edges on compact 'condensed' PAHs like circumcoronene, whereas 
trio-CH groups at 12.7 /jm arise at corners and thus imply more irregular PAHs. 

Heroic efforts in laboratory spectroscopy (e.g., Hudgins &: Allamandola 1999, see review by 
Tielens Sz Peeters 2004) have led to a large database of PAH spectra of various sizes and 
charges. It is clear that no single species fits all the features, but Allamandola et al. (1999) 
have succeeded in reproducing ISO spectra with a reasonable mixture of a dozen neutral and 
positively ionized PAHs, where each UIR feature is due to a superposition of many vibrational 
bands. The fact that the relative band strengths do not appear to vary much between PDRs 
exposed to orders of magnitude different radiation fields has been used to argue that the UIR 
features cannot be due to such an extended family of PAHs whose composition must change 
with physical conditions (e.g., Boulanger et al. 2000). However, as long as the ratio of UV 
intensity over electron density Go/n e (see §7) is similar and < 10 3 cm 3 , the spectra of a PAH 
family do not change much (Bakes et al. 2001, Li & Draine 2002). 

The high-quality spectra allow strict limits to be placed on functional groups attached to 
the PAHs, e.g., -CH3, -OH, -NH2 and -C=0 groups are at most 1% relative to aromatic C. 
The absence of strong 3.4, 6.85 and 7.3 /jm bands indicates that they do not contain large 
aliphatic moieties. The fact that the C-C stretch occurs at wavelengths as short as 6.2 ^m 
may be indicative of nitrogen substitution in the ring: pure PAHs start absorbing at 6.3 //m. 
Deuterated PAHs, in which an aromatic H is replaced by a D, have possibly been detected 
through bands at 4.4 and 4.65 /jm at high abundances, PAD/PAH> 0.1, in two sources (Peeters 
et al. 2004a). Other forms of large carbon molecules, e.g. Ceo or Cg~ , are not seen down to 
0.3% of total carbon (Moutou et al. 1999a). 

5.2. PAH Feature Variations 

Whereas the CH 3.3 and 11.2 fim bands correlate well with each other over a wide range of 
conditions, significant variations in the relative strengths of the C-H and C-C features are 
found, both between different types of sources and at different positions within a source (e.g., 
Hony et al. 2001). An excellent example of the latter case is provided by the M17 PDR, where 
the ratio of the 11.2/6.2 /im ratio changes by more than a factor of 2 across the ionization 
front (Verstraete et al. 1996). Since the C-H modes have much lower intrinsic band strengths 
for ionized PAHs compared with neutral PAHs, whereas the opposite holds for the C-C-modes 
(see Figure 9), this variation reflects a change in the degree of ionization: at the H II region 
interface, the PAHs are more ionized than deeper in the molecular cloud where they are also re- 
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Figure 9. Left: Source to source variations observed for the 6-9 /im PAH features 
with the SWS. Profile A is representative of all PDRs and reflection nebulae, including 
the diffuse interstellar medium. Profile B refers to the isolated Herbig Ae stars with 
disks, as well as most planetary nebulae. Profile C is observed for some post-AGB 
objects (Peeters et al. 2002a). Right: Laboratory spectra of a mixture of neutral 
PAHs compared to that of positively ionized PAHs. Note the stronger 3.3 //m feature 
for neutral PAHs and the stronger 6-9 /im complex for ionized PAHs (figure adapted 
from Allamandola et al. 1999). 

hydrogenated. The actual charge state depends on the ratio Go/n e and may include doubly- or 
even triply-ionized PAHs. The smaller PAHs may be destroyed by intense UV radiation inside 
the H II region (Verstraete et al. 2001), although observational evidence is still controversial 
(Abergel et al. 2003). 

The positions of individual C-H modes are relatively stable with peak wavelengths, changing 
by at most 4-6.5 cm -1 . In contrast, the C-C modes at 6-9 /jm vary by 25-50 cm -1 . Moreover, 
the variations in the different C-C profiles are linked to each other (see Figure 9, Peeters et al. 
2002a). Specifically, reflection nebulae and PDRs show a different 7.7 fim complex than isolated 
young stars and post-AGB stars. Some H II regions contain a 8.6 pLvn band which dwarfs the 
commonly much brighter 7.7 /xm band. While these changes clearly indicate variations in the 
composition of the PAH family in response to different physical conditions, the details remain 
to be understood, especially since few PAHs reveal strong spectral features longward of 7.7 /im 
in the laboratory. 

5.3. Hydrogenated Amorphous Carbon 

The C-H stretching mode of aliphatic (chain-like) hydrocarbons at 3.4 fim has been detected 
from the ground along many diffuse lines of sight with extinctions > 10 mag (e.g., Pendleton & 
Allamandola 2002). ISO has provided further contraints on the composition of this material by 
detecting new features at 6.85 and 7.25 /im toward the Galactic Center (Chiar et al. 2000). A 
good fit is obtained with Hydrogenated Amorphous Carbon (HAC) with H/C=0.5, with little 
evidence for strong C=0 modes which would be the relic of first-generation photoprocessed 
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ices. The total amount of carbon locked up in this material is uncertain, but may be as much 
as 25%. The relation with the aromatic material in the diffuse medium and the transformation 
between the two forms of carbon remains to be clarified (e.g., Mennella et al. 2001). 

6. SILICATES 

The detection of strong features at 9.7 and 18 fim bands due to the Si-0 stretch and O-Si-0 
bending modes (Gillett & Forrest 1973), coupled with the almost complete absence of silicon, 
magnesium and iron from the interstellar gas, convincingly demonstrates that a large fraction 
of the interstellar refractory material is in the form of silicates. ISO has allowed to (i) Obtain 
high quality silicate profiles of lines of sight through the diffuse interstellar medium; and (ii) 
Reveal a surprisingly rich mineralogy of crystalline silicates in material around young and old 
stars. 

6.1. Amorphous Silicates 

The 9.7 and 18 /im silicate bands are the deepest broad absorption features seen in ISO spectra 
toward embedded YSOs and background stars and are found to be very smooth. Demyk et al. 
(1999) fit the 9.7 and 18 /im bands for dense cloud material toward protostars with amorphous, 
porous pyroxenes (Mg^Fei-^SiOs) with some admixture of aluminum in the silicates. Kemper 
et al. (2004) obtain a good fit for the diffuse cloud 9.7 /j,m band toward Sgr A* with ~ 0.1 
/xm-sized amorphous silicates consisting mostly of amorphous olivine (~85%; MgFeSi04) with 
a small admixture of amorphous pyroxene (~15%). The limits on the degree of crystallinity 
are very strict, <l-2% in the case of dense clouds, and <0.4% for the more diffuse medium. 
Alternative fits by Bowey & Adamson (2002) with a complex mixture of crystalline material, 
constructed such that the many individual peaks blend into a smooth broad profile, seem less 
likely. 

This low degree of crystallinity is in stark contrast with the significant amounts of crystalline 
material of typically 10-20% (with extremes up to 75%) found in the envelopes around evolved 
stars (e.g., Molster et al. 2002). Since these stars are thought to provide a significant fraction 
of the interstellar dust, the observed lack of crystallinity has important implications. One 
explanation is that amorphization due to energetic particles occurs in the interstellar medium 
on a timescale of <10 Myr, significantly shorter than the dust destruction time scale. Another 
possibility is that the production rates of amorphous silicates from other sources (e.g., super- 
novae) are much larger than thought before, diluting the crystalline silicate fraction originating 
from (post-)AGB stars. 

6.2. Crystalline Silicates 

In dense clouds close to hot stars, the dust particles become sufficiently warm at a few hundred 
K to show silicate bands in emission. ISO examples are provided by Jones et al. (1999) for M 
17, Cesarsky et al. (2000a) in Orion, Lefloch et al. (2002) in the Trifid nebula and Onaka & 
Okada (2003) for Carina and S 171. Although the bulk of the emission is thought to be due to 
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Figure 10. Left: SWS spectra of disks around isolated Herbig Ae stars, illustrating 
their different spectral characteristics. HD 179218 is an example of a source with 
crystalline silicates, whereas HD 100453 does not show any 10 /xm emission. Both 
sources have rising SEDs (Group I). HD 144432 is an example of a source with a flat 
SED (Group II). The HD 144432 and HD 100453 spectra have been shifted by -5 
and —10 Jy for clarity. Right: Sketch of the different geometries for Group I and 
II sources. The silicate emission is thought to arise from the warm upper layers in 
the inner disk, whereas the PAH emission can also come from the flaring outer part 
(Meeus et al. 2001, Dullemond & Dominik 2004). 

small amorphous grains, some evidence for discrete features has been claimed. For example, 
Onaka & Okada find a feature at 65 /xm which may be due to diopside, a calcium-bearing 
crystalline silicate. Ceccarelli et al. (2002c) show that the spectral energy distribution of one 
protostar, NGC 1333 IRAS4, has excess emission at 95 /£m, which may be ascribed to calcite 
— a calcium-containing carbonate. Calcite has also been found in two planetary nebulae by 
Kemper et al. (2002) and is common in meteorites, with abundances of ~ 0.3 — 1% of the 
(warm) dust mass in all types of objects. If confirmed, the detection of calcite in a cold 
protostellar envelope would be additional evidence that its formation does not require liquid 
water. 

The richest silicate spectra are found in disks around young and old stars. A survey by 
Meeus et al. (2001) of a sample of isolated young intermediate mass stars away from molecular 
clouds — the so-called Herbig Ae stars (Waters & Waelkens 1998) — reveals silicate emission 
for at least 70% of the sources. In some objects, only the broad 10 and 18 fim bands due 
to amorphous silicates are seen, whereas others reveal a large variety of narrower solid-state 
emission bands (see Figure 10). The most spectacular example is provided by HD 100546, 
which shows at least 8 bands which can be ascribed to forsterite, Mg2Si04 (Waelkens et al. 
1996, Malfait et al. 1998). Other minerals identified in disks around Herbig Ae stars include 
enstatite (MgSiOa; Bouwman et al. 2001), hydrous silicates such as montmorillonite (Malfait 
et al. 1999), silica (Si02; Bouwman et al. 2001), metallic iron and iron oxide (van den Ancker 
et al. 2000a, Bouwman et al. 2000), although the FeO 23 /j,m feature can also be ascribed to 
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Figure 11 (color). Complete SWS-LWS spectrum of the isolated Herbig Ae star 
HD 142527, together with the decomposition into the different components, including 
hot and cold amorphous silicates, crystalline H2O ice, and tentatively the crystalline 
hydrous silicate montmorillonite (Malfait et al. 1999). 

FeS (Keller et al. 2002). The necessary laboratory spectroscopy is summarized in Henning 
(2003). Long wavelength data are often essential to make firm identifications (Figure 11), but 
LWS spectra exist for only a few sources. Typical mass fractions of the crystalline material 
are 5-10% compared with amorphous olivine, with a few % of silica. Several spectra also show 
features due to PAHs (see §5) and crystalline H2O ice (see §4), some of which may be blended 
with the crystalline silicate bands. Remarkably, no trend has been found between the age of 
the star and the presence of crystalline silicate or PAH material over the 0.1-10 Myr range. 

Although not related to age, the 10 (j,m. silicate profiles do show variations from source to 
source and are significantly different from those found in the diffuse medium (§6.1). In Figure 
12, a collection of silicate profiles is presented, arranged such that the peak emission occurs 
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Figure 12. Continuum-subtracted 10 /mi silicate emission bands observed toward a 
selected sample of isolated Herbig Ae stars with the SWS. For reference, the silicate 
bands of interstellar medium dust and comet C/1995 01 Hale-Bopp are included. The 
dashed line indicates the position of the 9.8 /tm amorphous silicate band observed in 
the interstellar medium. The full lines indicate the best-fit models using a mixture 
of amorphous and crystalline material with different grain sizes. The sharp 11.3 /im 
feature is evidence for crystallization; the shift toward longer wavelengths for grain 
growth (Bouwman et al. 2001). 



at progressively longer wavelengths. These shifts can be caused by two effects, collectively 
called 'dust processing': coagulation of grains resulting in larger average grain sizes, and crys- 
tallization resulting in discrete peaks, for example by forsterite at 11.3 p, Bouwman et al. 
(2001) argue on the basis of correlations between features that grain growth from sub-micron 
to micron-sized grains dominates the shift. The timescales for coagulation and annealing are 
not coupled, with crystallization having a longer timescale. The Si02 band at 8.6 /jm may be 
a by-product of the annealing process leading to crystalline forsterite, as has been observed in 
the laboratory (Rietmeijer et al. 2002). 

The absence of the 10 /im band in some sources is intriguing and may be due either to the 
absence of any small few /im) silicate grains or to low temperatures caused, for example, 
by shadowing (see §10). A recent re-analysis of the SWS spectrum of HD 100453 shows that 
both effects likely play a role (Meeus et al. 2002, Vandenbussche et al. 2004): features due to 
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large (~ 2 /im) crystalline forsterite grains are detected at longer (>30 /mi) wavelengths, but 
no emission is seen at 11 /im. 

Silicate emission from disks around low-mass T Tauri stars has been observed with PHOT-S 
(Natta et al. 2000, Giirtler et al. 1999). Although the quality of the spectra is limited, broad 
amorphous silicate features from (sub)-micron sized grains consisting of a mix of olivines and 
pyroxenes are found in all cases. There are hints of 11.3 /mi crystalline forsterite and 8.6 
lira, silica features in some sources, but confirmation requires higher S/N and higher spectral 
resolution data. Indeed, crystalline silicates in disks around solar-mass stars have recently 
been found from ground-based data by Honda et al. (2003). Mid-infrared spectroscopy of such 
objects is an obvious target for future space missions. 

7. PHOTON-DOMINATED REGIONS 

In Photon-Dominated or Photo-Dissociation Regions (PDRs), far-ultraviolet photons (6-13.6 
eV) control both the thermal and chemical structure of the neutral gas. Except for the densest 
shielded regions, the bulk of the molecular gas in our Galaxy and external galaxies has only 
a moderate opacity to UV radiation. Indeed, PDRs dominate the diffuse infrared emission 
from molecular clouds, through conversion of the ambient UV photons to infrared radiation 
from dust and gas. The study of PDRs was triggered by early airborne observations of strong 
far-infrared atomic fine-structure lines (e.g., Melnick et al. 1979). Prominent [C II] 158 /im and 
[O I] 63 and 145 /im lines from regions like Orion stimulated the development of detailed PDR 
models (Tielens k, Hollenbach 1985). Together with subsequent observations of UIR/PAH 
emission (see §5) and ground-based [C I] and high-excitation CO lines, a flurry of activity on 
PDRs resulted in the 1980's and 1990's. Comprehensive reviews have been given by Hollenbach 
& Tielens (1997, 1999). 

PDRs have been widely observed with ISO using all four instruments. CAM images at 7 /im 
(dominated by PAHs) and 15 /im (dominated by small grains) beautifully reveal the cloud 
surfaces set aglow by the surrounding UV radiation (see Fig. 1). Both CAM-CVF and PHOT- 
S have obtained spectra of the diffuse medium (Boulanger et al. 2000, Mattila et al. 1996). 
LWS and SWS spectra show a wealth of atomic fine-structure lines of neutral or singly ionized 
species, as well as the mid-infrared H 2 lines. Bright PDRs studied with ISO spectroscopically 
include S 140 (Timmermann et al. 1996), NGC 2023 (Moutou et al. 1999b, Draine & Bertoldi 
1999), NGC 7023 (Fuente et al. 1999, 2000), NGC 2024 (Giannini et al. 2000), W 49N (Vastel et 
al. 2001), S 106 (Schneider et al. 2003), S 125 (Aannestad & Emery 2001, 2003), S 171 (Okada 
et al. 2003), IC 63 (Thi et al. 1999), the Trifid nebula (Lefioch et al. 2002), and Trumpler 14 
(Brooks et al. 2004). 

Models have shown that the strength of the emission lines depends mainly on two parameters: 
the gas density n and the intensity of the incident radiation field, characterized by an enhance- 
ment factor Go compared with the standard interstellar radiation field. Here Go=l refers to 
an integrated ultraviolet (6-13.6 eV) intensity of 1.6 x 10 -3 erg s _1 cm -2 (Habing 1968), a 
factor of 1.7 lower than the reference field by Draine (1978). Most of the PDRs mentioned 
above have densities of at least 10 4 cm -3 and radiation fields Go > 10 3 . Since many of these 
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bright PDRs have also been observed with the KAO, the ISO data have helped mostly to refine 
the physical structure of the models (e.g., more accurate density determinations, constraints 
on clumpy or layered structures) and geometry (e.g., inclination, face-on versus edge-on, 3D 
geometry). Care should be taken not to overinterpret the data, however. Comparison of results 
from different PDR codes shows differences in predicted line intensities up to a factor of 2-3 
owing to different assumptions in the input data and chemistry. Also, some lines, in particular 
the [O I] 63 /im line, are optically thick so that model results are sensitive to the details of the 
radiative transfer and possible self-absorption by cold foreground gas. 

The higher spatial resolution ISO data have confirmed earlier suggestions that most of the 
observed [Si II] 34.8 pm emission originates from the outer PDR layer with Ay <2 mag rather 
than from any ionized gas in the region (e.g., Fuente et al. 2000, Lefloch et al. 2002), with the 
enhanced gas-phase silicon most likely due to photodesorption of Si-containing grain-mantle 
material (Walmsley et al. 1999). Not all PDRs show [Si II] emission, however (Rosenthal et 
al. 2000). 

The main new insights due to ISO have been (i) Extension of PDR studies to the lower n and Go 
regime; (ii) Direct measurement of the efficiency of heating the gas through the photoelectric 
effect; and (iii) Detailed studies of the H 2 mid- infrared pure rotational lines leading to new 
insights into the temperature structure and H 2 formation rate at high temperatures. PAHs 
are discussed mostly in §5. 

7.1. Low- Excitation PDRs and the Gas Heating Efficiency 

Low excitation PDRs with n < 10 4 cm -3 and/or G < 100 studied by ISO include L1457 
(MBM12) (Timmermann et al. 1998), the S 140 extended region (Li et al. 2002), L1721 (Habart 
et al. 2001), p Oph main (Liseau et al. 1999), p Oph W (Habart et al. 2003), and a set of 
translucent clouds also studied by optical absorption lines (Thi et al. 1999). Cederblad 201 
(Kemper et al. 1999, Cesarsky et al. 2000b) is an example of a PDR exposed to a cooler B9.5 
star where the radiation field has fewer far-UV photons. The physical structure of these clouds 
is often well characterized from complementary data, so that they form good tests of the basic 
PDR models. While some ISO lines (e.g., [C II]) are well explained, others (e.g., H2) cannot 
be reproduced with standard model parameters (see §7.2.3). 

Comparison of the total intensity of the main cooling lines (principly [C II] 158 ^m for the 
low-excitation PDRs) with the far-infrared continuum gives directly the heating efficiency e 
of the gas due to the photoelectric effect. Assuming that the lines are optically thin, values 
ranging from 1-3% are obtained. These efficiencies are somewhat higher than the range of 
0.1-1% found for a sample of higher n, Go PDRs by Stacey et al. (1991) with the KAO and 
that of 0.7% for S 106 studied with ISO. For very high values of Go/n appropriate for clouds 
like W 49N, the photoelectric heating efficiency is clearly suppressed by an order of magnitude 
to (0.01-0.04)%. This is consistent with models by Bakes & Tielens (1994), who show that 
intense radiation results in more highly charged grains which have a larger barrier for the 
photoelectrons to escape. So far, W 49N is the only PDR for which such low efficiencies have 
been found. Even Trumpler 14, a PDR powered by at least 13 O stars, and 30 Doradus in 
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Figure 13. Left: Pure rotational H 2 lines observed toward IC 63 with the SWS (Thi 
et al. 1999). Right: H2 excitation diagram for IC 63 compared with that derived from 
UV observations of the diffuse cloud ( Oph (Spitzer & Cochran 1973). The statistical 
weights are gj = 2 J + 1 with gi=l for para-H 2 and gi=3 for ortho-H 2 . 

the LMC, powered by at least 30 O stars, have efficiencies of 0.2-0.5% (Brooks et al. 2004). 
One possible explanation is that the ionization front in these cases is at such a large distance 
from the stellar cluster that the PDRs at the outer edges of the cloud do not have extreme 
parameters in terms of Go or n. 

Habart et al. (2001) went one step further by spatially correlating the cooling lines with 
infrared emission due to PAHs (4-12 A; IRAS 12 /an), very small grains (VSGs, 12-150 A; 
IRAS 25 and 60 /xm) and big grains (> 150 A; IRAS 100 //m). They conclusively show that the 
PAHs are the most efficient heating agents, with €pah=3%, evsG=l% an d €bg=0.1%. Habart 
et al. (2003) find £pah=4% for the p Oph West cloud. Both studies also find evidence for 
an increase in the PAH abundance near edge of cloud, consistent with other studies based on 
IRAS and/or ISO imaging (e.g., Boulanger et al. 1990, Abergel et al. 2002). 



7.2. H 2 Pure Rotational Lines 



7.2.1. Rotational temperatures. The H 2 pure rotational quadrupole lines are readily 
detected from bright PDRs, in spite of their small intrinsic Einstein A— values (see Figure 
13). In contrast with the vibration-rotation lines at 2 ^m seen from the ground, the lowest 
rotational levels are primarily excited by collisions in warm gas rather than by UV pumping 
(e.g., Black & van Dishoeck 1987, Sternberg & Dalgarno 1989, Draine & Bertoldi 1996). This 
is directly reflected in their excitation diagrams (see Figure 13, right), which show that the v=0 
rotational distribution is characterized by a temperature T rot =400-700 K (e.g., Timmermann 
et al. 1996, Fuente et al. 1999, Wright 2000), much lower than the vibrational temperatures 
of T vi b ~2000 K. The inferred T rot for various PDRs are remarkably similar in spite of their 
different characteristics. It should be noted, however, that only few sources have been studied 
in detail and that in several cases the integrations were not deep enough to reveal the lowest 
J = 2 — line at 28.22 fim. Where detected, this line shows evidence for a lower temperature 
component with T rot «100 K. It is interesting to note that the rotation diagrams derived 
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from H 2 UV absorption line data of diffuse (e.g., Spitzer & Cochran 1973, see Figure 13) 
and translucent (e.g., Snow et al. 2000) clouds are also similar. These data, which include 
the populations of the lowest two J=0 and J=l levels, clearly reveal the multi-temperature 
structure. Similar diagrams have been found for extragalactic regions in normal (e.g., Valentijn 
& van der Werf 1999) and starburst galaxies (e.g., Rigopoulou et al. 2002). 

In spite of their simple appearance, the interpretation of these diagrams has to take account of 
several factors. First, the gas kinetic temperature varies rapidly through the PDR layer from a 
value of several hundred K at the edge to less than 30 K when Ay > 1 mag. The pure rotational 
H2 lines arise primarily from the outer warm layer, but its temperature depends strongly on n 
and Go, as well as the assumed PAH heating parameters (e.g., Kaufman et al. 1999). Second, 
only the lowest rotational levels are collisionally excited. For diffuse clouds, this holds only up 
to J=2 whereas for denser PDRs like S 140 levels up to J=5 may be populated by collisions. 
Finally, the rotational populations are affected by UV pumping in at least two ways: in an 
absolute sense by providing more population in the higher-J levels, and in a relative sense by 
the fact that the excitation rates out of the various J levels differ with depth into the cloud. 
No PDR model has yet been published which fits all the ISO data and takes all these effects 
into account in a fully self-consistent manner. For example, the excellent fits to the S 140 H2 
data by Timmermann et al. (1996) use a fitted gas temperature profile rather than a computed 
one. 

One general conclusion of all PDR models is that the gas temperature needs to be higher in 
the zone where H 2 is formed than can be explained by the standard gas heating and cooling 
processes. Draine h Bertoldi (1999) discuss possible remedies. In addition to changes in 
geometry and density structure, possible uncertainties in the microphysics include the H 2 
inelastic collisional rate coefficients, the grain photoelectric heating efficiency, and an increased 
dust-to-gas ratio in the PDR zone due to anisotropic illumination of grains which causes them 
to drift through the gas. An additional option, namely an increased H 2 formation rate, is 
discussed below. Yet another explanation, namely that of weak shocks in diffuse gas, remains 
to be quantified. 

7.2.2. Ortho/para ratio. There has been significant confusion in the literature on the 
ortho/para ratio of H2, with several claims of ortho/para ratios that are lower than the equi- 
librium value corresponding to the temperature derived from excitation diagrams (e.g., Fuente 
et al. 1999). As pointed out by Sternberg & Neufeld (1999), care has to be taken with such 
statements since a low ortho/para ratio in excited levels naturally results from the fact that 
the optical depth in the UV pumping lines is higher for the ortho- than for the para-H2 lines 
(see their Figure 5). None of the ISO data measure the bulk of ortho- and para-H2, which 
resides in J=0 and 1. Additional concerns are mentioned by Wright (2000). Nevertheless, 
some of these claims of 'non-equilibrium' ortho/para ratios in high temperature gas may be 
valid (see also §8), in which case they are interesting indicators of dynamical processes, such 
as the advection of colder gas into the PDR. 

7.2.3. Increased H 2 formation rate. Detailed studies of several well-characterized, low- 
excitation PDRs have shown convincingly that the intensities of the lowest H 2 lines cannot be 
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reproduced with standard PDR parameters (e.g., Bertoldi 1997, Kemper et al. 1999, Thi et 
al. 1999, Habart et al. 2001, 2003, Li et al. 2002). The discrepancies are significant, at least 
an order of magnitude. Since for these regions the gas-grain drifts discussed above are not 
significant, the discussion has focussed on two solutions. First, an increase in the photoelectric 
heating rate, for example through a larger PAH or very small grain abundance. Second, an 
increase in the H2 formation rate at high temperatures. A combination of these two solutions 
has the effect of moving the H-H2 transition zone closer to the edge of the PDR where the 
temperatures are higher. The most quantitative study is by Habart et al. (2003a,b, 2004), 
who show that for most (but not all) sources an increase in the H 2 formation rate by a factor 
of ~ 5 at temperatures of a few hundred K is needed, compared with the standard rate of 
(3-4) xlO -17 cm 3 s -1 derived from UV absorption lines by Jura (1975) and Gry et al. (2002) 
for diffuse clouds. 

Since the residence times of physisorbed H on grains at such high temperatures are much too 
short, this conclusion has fundamental implications for our understanding of basic molecule 
formation processes. One possibility is that direct reactions of atomic H arriving from the gas 
on top of chemisorbed species dominate in an Eley-Rideal type process (see e.g., Duley 1996). 
Another option is that physisorbed H diffuses over the surface to find a chemisorbed H in a 
Langmuir-Hinshelwood-type process (Hollenbach & Salpeter 1971). Quantitative estimates of 
the rates of both processes have been made (Cazaux & Tielens 2002, Habart et al. 2004). The 
efficiencies depend sensitively on the number of chemisorbed sites and the activation energy 
barrier to molecule formation, both of which are uncertain. 

8. SHOCKS 

Shocks in molecular clouds can be caused by a large variety of processes including outflows and 
jets from young stars, supernovae and expanding H II regions. Recent reviews of the theory 
of shocks have been given by Draine & McKee (1993) and Hollenbach (1997). In this section, 
ISO observations of well-characterized shocks away from the driving sources are discussed; 
those in the inner parts of protostellar envelopes are included in §9. The continuum emission 
due to warm dust offset from the YSOs is generally very low so that all lines are in emission 
(see Figure 2). Early models by Hollenbach et al. (1989) show that slow, non-dissociative 
shocks have many strong diagnostic lines at infrared wavelengths. In these so-called C— type 
shocks, the hydro dynamical variables change continuously and the maximum temperature is 
only 2000-3000 K, too low to dissociate molecules. In contrast, the temperatures in high- 
velocity J— type shocks jump to more than 10 5 K, and most of the emission is from atomic 
lines at UV wavelengths. Also, these shocks can cause shattering or non-thermal sputtering of 
grain cores, resulting in enhanced abundances og Si, Mg, Fe and other grain constituents. 

ISO's main contributions to the study of shocks have been to (i) Better characterize the shock 
structure and physical conditions, especially that found in the slower shocks; (ii) Measure 
the total cooling power directly from observations, in particular the individual contributions 
from H2O and H2, and thus understand the overall energy budget; and (iii) Study global and 
systematic trends by observing much larger samples. A good overview has been given by 
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Figure 14. Left: narrow-band image of the L1448 outflow at 2.12 /im (H 2 +continuum), 
with superposed the LWS (long-dashed circles) and SWS (short-dished rectangles) aper- 
tures. The thick and thin solid lines delineate the contours of the blue- and red-shifted 
CO outflow from the L1448-mm source (Nisini et al. 2000). Right: LWS continuum- 
subtracted spectrum toward the Class source L1448-mm compared with the more 
evolved L1641N and LkHa 234 sources. Note the decrease of molecular emission, in 
particular H2O, toward the latter sources (Benedettini et al. 2003). 

Noriega-Crespo (2002). 

The ISO data have also convincingly shown that all broad band emission from shocks is due to 
lines (e.g., Cesarsky et al. 1999). For example, the ISOCAM 7 and 15 (im bands are dominated 
by the H 2 S(4) and S(l) lines, whereas the IRAS 12, 25, and 60 ^m emission is due to H 2 S(2), 
S(0), and [O I] 63 /im, respectively. 

8.1. Shock Structure and Physical Conditions 

Shocks associated with outflows studied by ISO include Orion (Harwit et al. 1998, Rosenthal et 
al. 2000, Sempere et al. 2000, Gonzalez- Alfonso et al. 2002), HH 1-2 (Cernicharo et al. 2000b, 
Lefloch et al. 2003, Cabrit et al. 2004), HH 7-11 (Molinari et al. 2000), HH 46 (Nisini et al. 
2002), HH 52-54 (Liseau et al. 1996, Nisini et al. 1996), HH 80-81 (Molinari et al. 2001), Cep 
A West (Wright et al. 1996, Froebrich et al. 2002), Cep E (Moro-Martm et al. 2001, Smith et 
al. 2003), DR 21 (Wright et al. 1997), L1448 (Nisini et al. 1999a, 2000; Froebrich et al. 2002), 
B 335 (Nisini et al. 1999b), L1551 (White et al. 2000), NGC 1333 (Bergin et al. 2003) and a 
variety of other low-mass YSOs including L1157 (Giannini et al. 2001). Shocks associated with 
supernova remnants investigated by ISO include IC 443 (Cesarsky et al. 1999) and 3C 391, W 
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Figure 15. Left: histogram of the total molecular cooling, defined as L mo i = Lh 2 o + 
Lqo + Lou, and of the cooling due to [O I], compared with the source bolometric 
luminosity for a set of embedded YSOs (Nisini et al. 2002). Right: Luminosity emitted 
by different molecular and atomic species along the L1448 outflow (see Figure 14 for 
positions). "Warm" and "hot" H2 refer to the pure rotational and vibration-rotation 
lines, respectively (Nisini et al. 2000) 

44 and W 28 (Rho & Reach 2003). In all types of shocks, molecular (H 2 , CO, H 2 0, OH) and 
atomic ([O I], [C II], [S I]) lines are often detected (see Figure 14 and 15), although there are 
sources in which no molecular lines except H 2 are seen (e.g., L1551 IRS5, HH 80-81). 

In spite of the fact that the available diagnostics often differ depending on whether SWS, 
LWS and/or CAM-CVF data are obtained, the conclusion is invariably that no single planar 
static shock can explain all the data. Specifically, the H2 excitation and the presence of other 
molecules require a combination of slow (<20 km s _1 ) and fast (30-50 km s _1 ) C— shocks with 
different covering factors, whereas the atomic lines can only be explained with a J— type shock 
(70 — 140 km s _1 ). The estimated filling factors of the shocks are small, typically a few % of the 
LWS beam. It is still unclear whether these different shocks can be combined in a single unified 
model. Attempts to do so include bow shocks with long flanks, where a broad distribution of 
shock speeds arises from the curvature (e.g., Froebrich et al. 2002), time-dependent C— shocks 
with embedded transient J— shocks (e.g., Chieze et al. 1998), or multiple shocks linked to 
multiple episodes of outflow activity (e.g., Nisini et al. 2000). 

8.1.1. Atomic lines, J-shock diagnostics. A clear diagnostic of a fast dissociative J- 
shock is the [Ne II] 12.8 /xm line, whose intensity is directly proportional to v s . This line is 
commonly associated with high-mass YSOs, but has been detected in only one low-mass YSO 
(HH 1-2: see Figure 2). The lack of [Ne II] or other highly ionized fine-structure lines toward 
low-mass YSOs, coupled with the detection of optical high-excitation lines, indicates that the 
area occupied by high velocity ~ 200 km s _1 shocks is small, typically < 1". Instead, van den 
Ancker et al. (2000b, c) show that [S I] 25.2 /im may be a better shock diagnostic. Another 
often- used probe is the [C II] /[O I] fine-structure line ratio, which is expected to be low in 
fast J-shocks. [C II] can also arise in PDRs, however, either from the surrounding molecular 
cloud emission or produced locally by the shock. Indeed, Molinari et al. (2001) and Molinari h, 
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Noriega-Crespo (2002) have argued that much of the observed [C II] emission in Herbig-Haro 
sources arises from PDRs at the walls of the outflow cavity. This radiation can also create a 
PDR in the quiescent gas ahead of the shock. 

8.1.2. H2 lines, C-shock diagnostics. As for PDRs, the pure rotational H2 lines dominate 
the mid-infrared spectra and are characterized by excitation temperatures of 700 — 1000 K. 
In contrast, the well-studied H 2 vibration-rotation lines usually give T ex =2000-3000 K. Many 
papers treat these two data sets separately, but the combined pure-rotation and vibration- 
rotation data observed with ISO clearly indicate that two excitation temperatures are common, 
with evidence for even higher excitation conditions (> 10000 K) in some sources (Cep A: Wright 
et al. 1996, Orion: Rosenthal et al. 2000). The most impressive collection of H 2 lines is seen 
for Orion at the so-called Peak 1 (Rosenthal et al. 2000, see Figure 3) and Peak 2 (CM Wright, 
private communication) positions, where 56 H2 lines have been seen including the v=0 S(l) 
to S(25) lines with upper energy levels up to 42500 K. The detection of the S(25) line is 
particularly intriguing and may be indicative of a non-thermal excitation processes, either by 
direct population of high- J levels during the formation of H 2 (so-called 'formation pumping') 
or by non-thermal collisions between ions and molecules in the shock. However, Le Bourlot 
et al. (2002) have recently re-evaluated the chemistry and excitation of H 2 in shocks using 
much improved H2-H and H2-H2 collisional rate coefficients, and concluded that the critical 
velocity at which H 2 is dissociated is significantly higher than previously thought. As a result, 
higher excitation of H 2 can result from thermal collisions only, and the Orion H 2 data can be 
reproduced with a two-component C— shock model without the need for non-thermal excitation 
processes. CAM-CVF images show that H2 avoids the high-excitation optical knots and has a 
distribution more consistent with shocked gas swept-up by the bow-shock (e.g., Cabrit et al. 
1999, 2004; Larsson et al. 2002) 

Theoretical calculations of the H 2 rotational excitation in various shock models give T ex values 
of 200-540 K for J— shocks, but a much larger range of 100-1500 K for C— shocks (see Figure 8 
in van den Ancker et al. 2000b). T ex does not depend much on density, but increases strongly 
with shock velocity. Thus, the observed range of T ex discussed above naturally leads to at least 
two C— type shocks. Interestingly, the rather narrow range of T ex predicted for J— shocks is 
similar to that found for PDRs (see §7.2.1). For PDRs, T ex does not depend much on density 
and is a good diagnostic of the UV radiation field Go- Except for the shock associated with 
T Tau (van den Ancker et al. 1999), however, no case of a fluorescent contribution to the H2 
shock excitation has been found. More recently, Cabrit et al. (2004) have stressed the use 
of the H 2 pure rotational lines to also constrain the transverse magnetic field strength and 
abundance of small grains, both of which affect the shock details. 

Since the collisional rate coefficients of H2 with H are significantly larger than those with 
H 2 or He, the H/H 2 ratio in the pre-shocked gas is critical in the models. In general, much 
better fits to the shock data are obtained if this ratio is ~0.1, at least two orders of magnitude 
larger than that expected in dark clouds. For most shocks, the inferred ortho/para ratio is 3, 
consistent with the warm conditions. The major exceptions are HH 54 (Neufeld et al. 1998) 
and HH 2 (Lefloch et al. 2003), where significantly lower ortho/para ratios of 1.2-1.6 are found. 
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Figure 16. Examples of high-excitation emission and absorption lines of H 2 for the 
Orion shock observed with the LWS (o-H 2 5o5 — 4 i4 at 99.4931 /im, p-H 2 5is — 4 i4 
at 95.6273 fim) and SWS (o-H 2 4 32 -3 03 at 40.6909 /tm, p-H 2 5 33 - 4 04 at 35.4716 
/im) Fabry- Perot spectrometers (Harwit et al. 1998, Wright et al. 2000). 

This low ratio is thought to be a relic of the earlier cold cloud stage, where the timescale of 
the shock passage has been too short to establish an equilibrium ortho/para ratio. Thus, the 
ortho/para ratio can be used as a chronometer, indicating for HH 54 that the gas has been 
warm for at most 5000 yr. Timmermann (1998) shows that the ortho/para ratio only attains 
the high-temperature equilibrium value of 3 for shock speeds greater than 20-25 km s _1 . 

8.1.3. CO, H 2 0, and OH lines, C-shock diagnostics: The LWS spectra of shocks 
associated with nearby low-mass YSOs often contain several CO lines of comparable strengths, 
ranging from J= 14-13 up to 25-24 (see Figure 1 and 14). The excitation temperature of CO 
is generally somewhat lower than that of H 2 . Comparison of the H 2 0/CO luminosity ratio 
with shock models requires slow (10-20 km s _1 ) C— type shocks. Since submillimeter CO line 
profiles often indicate the presence of much higher shock velocities, one interpretation is that 
the slow shocks are proceeding into a medium which has been put into motion by a previous 
generation of shocks. Indeed, Nisini et al. (1999, 2000) associate the strong molecular emission 
found for L1448-mm directly with the extremely-high velocity bullets seen in CO which move 
at speeds up to 200 km s _1 . 

The most spectacular high-mass example is again provided by Orion, where Harwit et al. (1998) 
detected many H 2 emission lines resolved with the Fabry- Perot (see Figure 16). The observed 
line strengths are well reproduced by the C— shock model of Kaufman & Neufeld (1996) in 
which most of the gas-phase oxygen is driven into H 2 leading to an abundance of ~ 5 x 10 -4 . 
At shorter wavelengths, the pure rotational H 2 lines go from emission into absorption, at 
least in the direction of Orion-IRc2 (Wright et al. 2000). Not all shocks associated with high- 
mass YSOs resemble Orion, however, even when scaled for the larger distances: H 2 emission 
lines are usually not detected with ISO, indicating much smaller filling factors (Boonman et 
al. 2003a, see also §9.2). 
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8.2. Gas Cooling in Shocks 



The complete infrared spectral coverage provided by the ISO spectrometers allows a direct 
determination of the contribution of each species to the cooling. Figure 15 summarizes the 
results for a survey of embedded low-mass YSOs by Nisini et al. (2000, 2002). In many 
cases, H2, CO and H2O have comparable contributions on-source with no single species clearly 
dominating. Off source, the H2O contribution may be less. Together, the molecular lines can 
contain up to an order of magnitude more flux than the [O I] 63 //m line in the densest regions. 
The total observed cooling is often more than 50% of the mechanical power input of the shock, 
illustrating that the bulk of the kinetic energy is transformed into infrared radiation. 

In J— type shocks such as found in Herbig-Haro objects, the [O I] 63 ^m line is usually the 
dominant coolant. In that case, its intensity is a direct measure of the mechanical luminosity, 
which in turn is related to the mass-loss rate. Values of 10 -7 — 10 -4 M Q yr _1 have been 
inferred this way, which correlate well with those derived from CO outflow maps indicating 
that the two phenomena are physically associated (Saraceno et al. 1999, Liseau et al. 1997). 

9. EMBEDDED YOUNG STELLAR OBJECTS 

ISO spectra taken at the position of the protostars themselves often show most of the spectral 
characteristics of cold and warm gas and dust discussed above, arising from different parts 
of the circumstellar environment (see Figure 6, van Dishoeck et al. 1999). The cold outer 
envelope is seen prominently in ice absorption, whereas the PDRs and shocks are revealed 
through atomic, molecular or PAH emission lines. As the protostellar evolution progresses, 
the relative contribution of each of these components changes. The main advances made 
possible by ISO have been to (i) Quantify the relative contributions from different physical 
components; and (ii) Establish a protostellar classification scheme based on mid- and far- 
infrared spectroscopic features. For low- mass YSOs, this scheme complements and extends 
the Class 0, I and II classification based on the spectral energy distribution (SED) (e.g., Lada 
1999). For higher mass objects where the temporal ordering is much less clear (Churchwell 
2002), infrared spectroscopy is a powerful tool to characterize the earliest stages. 

For sensitivity reasons, low-mass YSOs have been studied primarily with the LWS and CAM- 
CVF, whereas high-mass YSOs have been intensively observed with the SWS. These two sets 
of sources are therefore discussed separately, but within a common physical and chemical 
framework. 

9.1. Low-Mass YSOs 

9.1.1. LWS surveys. Giannini et al. (2001), Nisini et al. (2002) and Benedettini et al. 
(2003) summarize the LWS observations of 28 low luminosity (<100 L ) YSOs, consisting of 
17 Class and 11 Class I objects. Several of these objects have been studied individually in 
earlier papers, e.g., IRAS 16293 -2422 (Ceccarelli et al. 1998a), B 335 (Nisini et al. 1999b), 
L1448-mm (Nisini et al. 2000), L1551 IRS5 (White et al. 2000), Serpens SMM1 (Larsson et 
al. 2002) and NGC 1333 IRAS 4 (Maret et al. 2002). The atomic [O I] 63 /jm and [C II] 158 
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jum lines are detected in all sources and molecular lines in many cases, most prominently in 
the Class objects. Specifically, CO, H 2 and OH are detected in 94%, 70% and 41% of the 
Class sources, respectively, whereas CO and OH are seen in 64% and 27% Class I objects. 
No H 2 has been detected in any Class I source, except for T Tau (see below). The spectral 
evolutionary sequence is illustrated in Figure 14. 

The total molecular cooling in Class objects is significantly larger than that in Class I sources, 
whereas the atomic cooling through [O I] 63 /ira is similar (see Figure 15). The total atomic + 
molecular line cooling Li; ne /Lboi ~ 10 -2 for Class sources, whereas it is an order of magnitude 
lower for Class I objects. This decline is consistent with that in the outflow momentum flux, 
which in turn is thought to be related to the decay of the mass accretion activity with time 
(Andre et al. 2000). 

Temperatures, densities, and column densities have been obtained by fitting molecular excita- 
tion models to the CO line ratios, giving values of a few hundred to >1000 K and densities 
of 10 5 — 10 6 cm -3 . The filling factor of this warm and dense gas is small, ~ 10 -9 sr or a few 
hundred AU radius. Assuming that the physical parameters derived from CO are also valid for 
H2O and OH, the inferred H 2 abundances range from 10 -5 to a few xlO -4 and increase with 
the gas temperature. Where detected, the OH abundance is at least an order of magnitide 
lower than that of H 2 0, of order 10 -6 — 10 -5 . For Class I objects, the upper limit on the 
H 2 abundance is < 10 -5 , indicating that the absence of H 2 is not simply due to a lower 
envelope mass, but is likely caused by photodissociation of H 2 to OH in the more diffuse, 
lower density Class I envelopes. Also, the longer timescale can result in H 2 freeze-out onto 
grain mantles after ~ 10 5 yr (Bergin et al. 1998). 

Two scenarios have been put forward for the origin of the molecular emission from Class 
objects. Ceccarelli et al. (1999, 2000) and Maret et al. (2002) favor a model of a collapsing 
envelope in which the lower-lying lines originate from the outer envelope and the higher-lying 
lines from the inner 'hot core' where ices evaporate (see Figure 6). In contrast, Giannini et al. 
(2001), Benedettini et al. (2002) and Nisini et al. (2000, 2002) argue that most of the emission 
arises in the J— and C— shocks associated with the outflow impacting the envelope, much 
like discussed in §8. Arguments in favor of the latter explanation include (i) the fact that in 
some cases the lines have similar strengths at the source and off-source outflow positions; (ii) 
the broad ~20 km s _1 H 2 line profiles seen by SWAS; and (iii) the linear relation between 
total line luminosity and the kinetic luminosity of the outflow derived from mapping of CO 
millimeter lines. The correlation with other molecular shock tracers such as SiO is less clear. 

It should be kept in mind, however, that the observed amount of warm gas is typically only 
0.01 Mq (a few % of the total envelope mass) and that evidence for warm inner envelopes 
and ice evaporation is clearly established from other data (see van Dishoeck 2003 for review). 
Larsson et al. (2002) argue that such models cannot reproduce the CO emission for one source, 
but the case for H 2 is less clear. For example, Maret et al. (2002) show that for NGC 1333 
IRAS4 — a source where the molecular emission peaks at the YSO — , the H 2 lines can be 
well explained with an abundance of 5 x 10 -7 in the outer region and 5 x 10 -6 in the inner 
envelope. For high-mass YSOs, ice evaporation is also known to contribute to the SWAS and 
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ISO H 2 observations (see §3.4 and 9.2). Finally, Ceccarelli et al. (2002c) suggest that some of 
the high-excitation CO may originate from a super-heated disk in the embedded phase. Future 
high spectral and spatial data with the Herschel Space Observatory are needed to determine 
the relative contributions from each of these scenarios for a large sample of sources. 

9.1.2. CAM-CVF and SWS surveys. As noted in §4.2, CAM-CVF observations, together 
with a few SWS spectra, reveal the presence of abundant ices in the cold outer envelope of 
low-mass YSOs. However, the lack of spectral resolution prevents the establishment of any 
trends based on changing line profiles or gas/solid ratios, such as observed for high-mass YSOs. 
Medium resolution mid-infrared spectroscopy of low-mass YSOs is an area that is wide open 
for future studies. 

9.1.3. T Tau. A special case is formed by T Tau, the prototypical Class II object but whose 
infrared line spectra are more characteristic of an embedded Class I object. Indeed, T Tau has 
one of the richest SWS (van den Ancker et al. 1999) and LWS (Spinoglio et al. 2000) spectra, 
full of atomic and molecular (H2, CO, H2O and OH) lines. In addition to shocks, the influence 
of UV radiation is clearly seen in the fluorescent H2 component and in the H2O-OH chemistry. 
T Tau may be peculiar because it is a close binary system with a dusty envelope seen face-on, 
filling a large fraction of the beam. 

9.1.4. FU Orionis objects: Another interesting sample is formed by the FU Orionis objects, 
studied with the LWS by Lorenzetti et al. (2000). These are low luminosity young stars which 
show a sudden increase in brightness, thought to be caused by a large increase in the disk 
accretion rate. In addition to the usual atomic [O I] and [C II] fine-structure lines, these 
objects surprisingly show [N II] 122 /im emission, not seen toward any other type of YSO. [N 
II] must arise from ionized gas, and since no higher ionization stages of other elements are seen, 
this gas must have a low ionization fraction and low density (n(e) < 100 cm -3 ). One possibility 
is that it is produced by enhanced UV radiation from the disk-star boundary layer. 

9.2. Intermediate and High-Mass YSOs 

Several samples of intermediate and high-mass sources have been studied with ISO, for different 
purposes. Lorenzetti et al. (2002) summarize the LWS observations on Herbig Ae/Be stars, 
covering luminosities from 5 to 10 5 L Q (see also Giannini et al. 1999, Lorenzetti et al. 1999). 
Most of these objects are associated with diffuse nebulosity and show [C II] and [O I] fine- 
structure emission characteristic of PDRs. Some of the extended PAH emission has been 
imaged by Siebenmorgen et al. (2000). High-luminosity (> 10 4 L Q ) deeply embedded sources 
studied with the SWS are summarized by Gibb et al. (2004), whereas a catalogue of SWS- 
LWS spectra of H II regions is presented by Peeters et al. (2002b). LWS spectra of other 
high-luminosity objects are summarized by Benedettini et al. (2003). A common characteristic 
of all these luminous sources is that they show almost no molecular emission lines, and often 
just a handful of atomic emission lines, usually only [O I] 63 /im and [C II] 158 /im arising 
in a PDR component. If an H II region has started to develop, higher ionization stages are 
seen as well and can be used to determine the effective temperature and spectral type of the 
embedded source. Molecular material, both in the gas and solid state, is however present and 
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Figure 17. SWS spectra of young stars at different stages of evolution (van den 
Ancker et al. 2000a,c; Malfait et al. 1998). From top to bottom — in a rough 
evolutionary sequence — the spectra change from dominated by solid-state absorption 
features and gas (shock) emission lines, to PAH features and PDR lines, to amorphous 
and crystalline silicates. 



is prominently seen in SWS absorption spectra. 

9.2.1. From shocks to PDRs. Two particularly illustrative examples of the combined 
SWS and LWS spectral evolution have been provided by van den Ancker et al. (2000b, c). The 
first concerns the BD +40°4124 group containing three YSOs with infrared luminosities of 
270-1600 L Q (see Figure 17). The LkHa 225 spectrum is characteristic of a deeply embedded 
source, with deep silicate and ice absorptions. LkHa 224 is an example of a new class of nearly 
featureless sources, showing no silicate or ice absorption or emission. The spectrum of BD 
+40°4124 itself, a B2Ve star, has strong PAHs and many atomic lines. Van den Ancker et al. 
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(2000c) show that the atomic and H2 lines for the embedded sources originate from shocks, 
whereas those for BD +40°4124 are characteristic of a PDR, with a contribution from the H 
II region. Thus, as the envelope disperses, the PDR component starts to dominate compared 
with the shock emission. 

A similar picture is found from comparing two high luminosity (2 — 4 x 10 4 L ) objects, Cep 
A and S 106 (van den Ancker et al. 2000b). The Cep A spectrum shows deep ice and silicate 
absorptions, but no PAH emission nor ionized atomic lines, indicating that it has not yet 
started to develop its H II region. In contrast, the S 106 spectrum is dominated by PAHs 
and by atomic lines from a wide range of ionization stages characteristic of a young 08 star. 
The atomic and molecular H2 and CO emission lines indicate shock excitation for Cep A and 
PDR excitation for S 106. The situation for other high-mass YSOs is often less clear-cut. For 
example, the Orion SWS spectra reveal a mixture of many physical components in the beam, 
including the foreground H II region, the PDR, the powerful outflow and the embedded YSO 
(van Dishoeck et al. 1998). 

9.2.2. Evolution within the deeply embedded stage. While the above characteristics of 
the deeply embedded versus the more evolved PDR/H II region stages are readily recognized, 
more subtle evolutionary effects within the deeply embedded stage can be derived from careful 
analysis of the ice and gas-phase absorption data toward high-mass protostars. These are 
very young objects, before the onset of the ultracompact H II region or 'hot core' stage. As 
discussed in §4.3, the profiles of several ice bands show strong evidence of progressive heating 
of the envelope. Similarly, the gas/solid ratios of species like H 2 increase systematically for 
the same sources, as does the color temperature of the cold dust. The most likely scenario, 
put forward by van der Tak et al. (2000b), is that the warmer sources have a lower ratio of 
envelope mass M env compared with the bolometric luminosity of the source Lboi, so that the 
overall temperatures are lower. If the lower M env are due to the gradual dispersion of the 
envelope rather than to initial conditions, the sources can be put in a time sequence with the 
warmer sources being more evolved, thus providing a powerful diagnostic of the first few xlO 4 
yr of protostellar evolution. 

Why do high-mass sources not show more prominent molecular emission lines in the SWS and 
LWS spectra, like their low-mass counterparts? The main reason must be beam dilution, since 
the high-mass sources often have ten times larger distances. Detailed modeling shows that the 
'hot core' region where the temperature is above that for ice evaporation of 90 K is typically 
1000 AU, less than 1" at >1 kpc. Apparently, shocks have an equally low filling factor, as 
evidenced from the absence of H2 or [S I] emission lines. Moreover, the far-infrared continuum 
of high-mass YSOs is so strong that weak emission lines are more difficult to detect at low 
spectral resolution. 

The case of H 2 is instructive of the abundance and excitation variations within high- mass 
envelopes. Boonman et al. (2003a) analyzed all available H 2 SWS, LWS and SWAS data 
in terms of possible abundance profiles. The SWS data indicate ice evaporation in the inner 
'hot core' with abundances of > 10~ 4 , whereas the absence of LWS lines and the weak SWAS 
emission require significant H 2 freeze-out in the outer cold envelope with abundances < 10~ 8 
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at T < 90 K (see also Snell et al. 2000). The SWAS profiles consist of narrow and broad 
components, with the outflow contribution at most 50%. In some cases, only a narrow (<few 
km s _1 ) profile is seen with a line strength consistent with optically thick H 2 emission at 
~ 100 K arising from a ~ 1" source. This emphasizes the need for future high spatial and 
spectral resolution H 2 observations to properly disentangle the shock, envelope and 'hot core' 
components. 

10. CIRCUMSTELLAR DISKS AROUND YOUNG STARS 

10.1. Dust Evolution and Disk Structure 

As discussed in §5 and 6, the mid- and far-infrared spectra of disks around young stars exhibit a 
wide variety of features due to crystalline silicates, PAHs and ices (Meeus et al. 2001, Bouwman 
et al. 2001). Although no trend with age within the pre-main sequence stage has been found, 
there are clear indications that the dust evolves from the diffuse cloud to the disk phase. 
Analysis of the 10 /jm silicate band shows two effects. The first is coagulation of 0.1 /im to 
micron-sized grains. Interestingly, small PAHs remain present during this process, although 
their spectral characteristics are clearly different from those of molecular clouds (Peeters et 
al. 2002a). The second effect is crystallization, resulting in an increasing fraction of emission 
from forsterite and other minerals when going from the Herbig Ae stars to a debris disk like 
(3 Pictoris and comets like Halley and Hale-Bopp (Figure 12). The near-complete absence 
of crystalline silicates in interstellar clouds indicates that the crystallization must occur in 
the disks themselves. The correlation of forsterite with silica suggests that thermal annealing 
is the dominant process, but this has to occur at temperatures greater than ~1100 K for the 
timescale to be shorter than the disk lifetime. Since the crystalline material is observed at lower 
temperatures further out in the disk, this implies that some radial mixing process must occur 
(see also Bockelee-Morvan et al. 2002), although the mechanisms are not yet well understood. 
The data indicate that the timescale for annealing is longer than that of coagulation. The 
origin of crystalline water in cold disks like that around HD 142527 (Figure 11) is also unclear. 

The mid infrared spectra, combined with their SEDs, also provide insight into the physical 
structure of disks. Natta et al. (2000) noted that the silicate emission seen in the PHOT-S 
spectra of T Tauri stars requires the presence of a disk 'atmosphere', in which the disk surface 
is heated from the top by radiation from the star (Calvet et al. 1992, Chiang & Goldreich 1997, 
Chiang et al. 2001). Meeus et al. (2001) divided the isolated Herbig Ae star ISO spectra into 
two groups: those with a rising SED (Group I) and those with a flat continuum (Group II). 
The SEDs of the latter group can be fitted with a power-law provided by an optically thick 
but geometrically thin disk, whereas the former requires addition of a black-body component 
with 100-200 K provided by a disk atmosphere in the inner ~ 10 AU (see Figure 10). 

Natta et al. (2001) and Dullemond et al. (2001) (see also Dominik et al. 2003) propose a 
slight modification, in which the inner disk is truncated at the dust sublimation temperature 
of ~1400 K at ~0.5 AU, resulting in a sharp rim exposed to intense radiation and producing 
most of the near-infrared radiation. Due to the higher temperatures, the rim will have an 
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increased scale height ('puffed up'), shielding the outer parts of the disk and resulting in 
lower temperatures in shadowed regions. Group I sources have a flaring outer disk directly 
exposed to the stellar radiation, where the silicate and PAH emission is thought to arise. In 
Group II sources, the disk may be so small that there is no flaring region. Alternatively, 
grain growth and/or dust settling may cause the disk to become self-shadowing (Dullemond & 
Dominik 2004) . The precise origin of the silicate emission in Group II sources is not yet fully 
understood. The presence of an additional spherical (remnant) envelope cannot be excluded 
for some cases (e.g., Miroshnichenko et al. 1999), and spatially resolved data will be essential 
to disentangle the different components. 

The disk around HD 100546 is clearly different from that around other stars: nearly 70% 
of the absorbed and re-emitted starlight of HD 100546 emerges at mid- rather than near- 
infrared wavelengths, indicating that the disk must intercept the radiation at larger distances, 
around 10 AU. Another peculiarity of this object is that the crystalline silicates have a lower 
temperature than the amorphous material and an increasing abundance with radial distance. 
One possibility is that a proto-Jovian planet has cleared out a gap around 10 AU, leading to 
a 'puffed-up' rim on the far-side of the gap (Bouwman et al. 2003). The wealth of crystalline 
silicates could have been produced in a collisional cascade of asteroidal-sized objects stirred 
up by the proto-Jupiter, or by shocks caused by tidal interactions of the planet with the disk 
producing crystalline material through flash heating. Some of this material could have been 
lifted and scattered to the disk surface and driven out by radiation pressure to larger distances 
High-spatial resolution data are needed to confirm this scenario, but it illustrates the potential 
of mid-infrared spectroscopy to trace planet-formation processes. 

The similarity of the HD 100546 and comet Hale-Bopp spectra (Crovisier et al. 1997) has 
often been cited as evidence of a strong interstellar medium-solar system connection, although 
the scenario described above is quite specific for one object. Additional evidence for such a 
connection has been derived from the fact that the ice abundances found in protostellar regions 
(see §4.2) are comparable to those observed in cometary comae (e.g., Ehrenfreund et al. 1997b, 
Bockelee-Morvan et al. 2000). This ice similarity is strong evidence that some fraction of the 
interstellar material must have been included largely unaltered, but the origin of the crystalline 
silicates in comets is still a mystery in this context. 

10.2. Gas in Disks 

An equally important issue is the timescale over which the gas disappears from the disk, 
since this directly affects the formation of gaseous giant planets. Ground-based millimeter CO 
surveys are hampered by freeze-out in the cold parts of the disk and photodissociation in the 
upper layers (e.g., Aikawa et al. 2002), whereas infrared CO lines probe only a small fraction 
of warm gas (e.g., Najita et al. 2003). The lowest H 2 pure-rotational lines provide a powerful 
alternative gas probe because it is the dominant constituent of the gas, does not freeze-out 
onto grains, self-shields effectively against photodissociation, and has optically thin lines. The 
main disadvantage is that the pure rotational lines are intrinsically weak and difficult to detect 
superposed on a strong continuum at low spectral resolution. Also, the gas temperature needs 
to be sufficiently high, >80 K. 
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Deep searches for the H2 S(0) and S(l) lines at 28 and 17 /jm have been performed with 
the SWS toward a dozen young stars (Thi et al. 2001a, b). These sources were chosen to be 
isolated objects to minimize contamination from surrounding cloud material. Detections of 
lines have been claimed in several sources, including potentially some older Vega-type objects, 
but confirmation with higher spatial and spectral resolution data is needed (see Richter et al. 
2002). Since the mid-infrared H2 lines are sensitive to a fraction of Mj up of warm gas, they are 
well worth pursuing with future instrumentation. Other gas tracers such as the [C II] 157.7 
/j,m line have been considered as well (Kamp et al. 2003). 

11. CONCLUDING REMARKS 

ISO has clearly demonstrated the power of full mid- and far-infrared spectra to probe a wide 
variety of processes associated with star- and planet formation, providing insight that cannot 
be obtained from imaging data alone. Using ISO, the inventory of ices has been completed, a 
census of the PAHs features has been obtained, a new diagnostic of warm gas -H2- has been 
opened up, key gas-phase molecules (H2O, CO2) have been widely observed, and a wealth of 
crystalline silicate features has been revealed. More importantly, ISO has shown systematic 
variations in the features during the star-formation process, providing a spectral classification 
scheme of the earliest embedded stages which are difficult to probe otherwise and stimulating 
new models. Related variations in the composition of the material when it is cycled from 
one phase to another are also found: freeze-out and evaporation of ices, ionization and re- 
hydrogenation of PAHs, and amorphization, crystallization and coagulation of silicates. At all 
stages, specific lines and line ratios provide powerful diagnostics of temperatures and densities, 
as well as the relative importance of shock- and photon heating of the gas. Since many of these 
same features and processes are also ubiquitous in extragalactic sources out to the highest 
redshifts, a good understanding of their diagnostic capabilities locally is essential. The ISO 
archive will remain a valuable — and in many ways unique — resource for future studies. 

At the time of writing, the first results from the Spitzer Space Telescope are just emerging. 
With its increased sensitivity up to two orders of magnitude, Spitzer will extend mid- and far- 
infrared imaging and spectroscopy into new regimes that ISO could not probe, such as medium 
resolution spectroscopy of embedded low-mass protostars, of disks around young stars down 
to the brown-dwarf regime, of debris disks, and of H2 and other lines in the most diffuse PDRs 
and weakest shocks. The ISO spectra summarized here form a guideline to the interpretation 
of the Spitzer data. Spitzer's spectroscopic capabilities are limited, however, both in spectral 
resolution and wavelength coverage. A combination of several future missions, in particular 
the Stratospheric Observatory for Infrared Astronomy (SOFIA), ASTRO-F, the Herschel Space 
Observatory, and the James Webb Space Telescope (JWST), will be needed to obtain spectra 
over the full 2.4-200 /im range covered by ISO. The highest spectral and spatial resolution 
data will be provided at mid-infrared wavelengths by large optical ground-based telescopes, 
and at submillimeter wavelengths by the Atacama Large Millimeter Array (ALMA) and other 
millimeter arrays. ISO has paved the way for the exciting science to be performed with these 
new facilities. 
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